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Preface 


Raman  scattering  studies  were  performed  on  III-V  Nitride  materials  including  binary 
crystals,  ternary  alloys,  and  superlattices  of  GaN/AIN  and  GaN/AlGaN.  Our  results 
show  that  in  AlGaN  alloys  the  asymmetric  behavior  of  the  E21ineshape  arises  from 
random  disorder;  the  E2  phonon  does  not  follow  a  one  mode  behavior  and  the  AlGaN 
films  on  SiC  are  under  small  tensile  stress.  In  the  InGaN  too  the  A1  mode  shows  a  one 
mode  behavior  while  two  mode  behavior  may  be  possible  for  E2.  Possible  mode 
mixing  takes  place.  Evidence  of  inhomogeneities  and  spinodal  decomposition  is  seen  in 
these  films.  Phonon  lifetime  measurements  are  made  in  the  binary  crystals  and  in  every 
case  we  see  evidence  of  two  lifetime  regions.  In  the  superlattices  we  see  anomalous 
low  temperature  behavior  that  is  consistent  with  a  density  of  states  of  spectral 
distribution  that  is  due  to  the  roughness  of  the  interfaces  and  which  agrees  with  TEM 
measurements.  Finally  preliminary  measurements  demonstrate  strong  interface 
phonons  in  the  GaN/AIN  superlattices  which  give  good  agreement  with  theory. 
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ABSTRACT 

A  detailed  study  of  photoluminescence  (PL)  of  GaN(lnm)/Al2oGa8oN(3.3nm)  twenty 
periods  superlattice  grown  via  metalorganic  chemical  vapor  deposition  is  presented.  The 
dependence  of  the  PL  emission  energy,  linewidth,  and  intensity  on  temperature,  in  the  low 
temperature  regime,  is  consistent  with  recombination  mechanisms  involving  bandtail  states 
attributed  to  a  small  degree  of  interfacial  disorder.  The  activation  energy  of  the  non-radiative 
centers  in  our  superlattice  agrees  well  with  the  value  we  derive  for  the  width  of  the  tail-state 
distribution.  Moreover,  we  find  that  the  average  phonon  energy  of  the  phonons  that  control  the 
interband  PL  energy  at  high  temperatures  is  larger  for  the  superlattice  than  for  a  high  quality 
GaN  film.  This  observation  is  consistent  with  model  calculations  predicting  the  phonon  mode 
properties  of  GaN-AIN  -  based  wurtzite  heterostructures. 
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Anomalous  temperature  dependence  of  the  photoluminescence  (PL)  emission  energy 
from  InGaN  quantum  wells  that  does  not  follow  at  low  temperatures,  the  accepted  Varshni  [1]  or 
Bose-Einstein  [2-3]  type  behavior  of  a  semiconductor  band-gap  has  been  reported  [4].  Such 
behavior  has  been  observed  as  well  in  quantum  wells  and  superlattices  of  the  zincblende 
structure  [5-8].  The  PL  energy,  instead  of  being  approximately  constant  at  low  temperatures, 
exhibits  a  red  shift  follow  by  a  blue  shift.  In  the  above  studies  the  anomalous  dependence  was 
attributed  to  the  effect  of  the  interface  defect  density  of  states  (DOS),  on  the  PL  energy. 

The  presence  of  defect  density  in  an  heterostructure  device  is  crucial  to  the  device 
performance,  especially  in  structures  where  multiple  interfaces  exist.  In  this  letter,  we  present  a 
detailed  study  of  the  photoluminescence  from  a  GaN/Alo^Gao.sN  superlattice  of  wurtzite 
structure.  We  establish  the  presence  of  a  defect  density  of  states,  which  controls  the  PL  energy  at 
low  temperatures  and  acts  as  a  non-radiative  channel  to  the  PL  at  higher  temperatures.  Moreover, 
the  electron-phonon  interaction  that  affects  the  PL  at  higher  temperatures  is  stronger  in  the 
superlattice  than  in  a  GaN  film  indicating  the  contributions  of  phonons  of  higher  energies. 

The  luminescence  measurements  were  acquired  utilizing  the  244  nm  (5.08  eV)  line  from 
the  second  harmonic  generation  of  the  488  nm  line  of  Ar  ion  Coherent  laser  and  the  Triplemate 
ISA  monochromator  in  conjunction  with  a  liquid  nitrogen-cooled  CCD  detector.  The  cold  stage 
consists  of  a  closed-cycle  cryostat,  which  enables  measurements  from  10K  up  to  room 
temperature.  The  superlattice  consists  of  twenty  periods  of  (lnm)GaN/(3.3nm)Alo.2Gao.gN  grown 
on  2pm  GaN  with  a  25nm  GaN  buffer  layer  on  a  (0001)  sapphire  substrate.  The  nitride  materials 
described  in  this  work  are  grown  by  metalorganic  chemical  vapor  deposition  (MOCVD)  in  an 
EMCORE  D125  UTM  rotating-disk  reactor.  The  growth  temperature  is  -1050  °C  for  the  thick 
undoped  GaN  layer  and  the  Alo.2Gao.8N  superlattice  layers. 

In  the  following  we  present  our  study  focusing  on  the  photoluminescence  behavior  as  a 
function  of  temperature  of  the  GaN/Alo^Gao.gN  superlattice.  Figure  1  depicts  the  PL  spectrum  of 
the  SL  acquired  at  10K.  The  peak  at  3.795  eV  with  a  full  width  at  half  maximum  of  32  meV  is 
attributed  to  the  PL  emission  from  the  GaN  wells.  The  additional  PL  peaks  observed  in  the 
spectrum  are  at  ~  3.488  eV,  due  to  exitonic  emission  of  the  underlying  GaN  buffer  layer,  as  well 
as  a  weak  peak  at  ~  3.7  eV,  the  origin  of  which  will  be  addressed  in  a  later  study.  The  inset  to 
Figure  1  depicts  the  intensity  of  the  luminescence  from  the  superlattice  as  a  function  of 
excitation  intensity  at  10  K.  A  linear  relation  between  the  two  intensities  was  ascertained  which 


is  indicative  that  the  PL  emission  of  the  SL  at  10  K  is  due  to  excitonic  recombination  [9].  Figure 
2  shows  the  luminescence  peak  position  as  a  function  of  temperature  of  the  superlattice 
compared  to  that  of  a  high  quality  GaN  thin  film.  The  inset  in  fig.  2  plots  the  PL  linewidths  of 
from  the  superlattice  sample. 

It  has  been  established  that  the  temperature  dependence  of  an  interband  transition  energy, 
E(T),  in  semiconductors  can  be  approximated  via  the  relation  [2-3,  10-12]: 

E{T)  =  E(o)  -  (1) 

This  relation  represents  the  modification  of  a  band-gap  of  a  semiconductor  due  to  the  electron- 
phonon  interaction  at  thermal  equilibrium,  and  may  be  derived  from  the  Bose-Einstein 
approximation  for  the  lattice  vibration  energy:  (E)  =  ( 1  / 2  +  («) )hco  where(/i)is  the  average 

phonon  occupation  number  [12],  In  Equation  1,  E(0)  is  the  transition  energy  at  0  K,  a  is  a 
measure  of  the  strength  of  the  electron-average  phonon  interaction,  and  ©  corresponds  to  a  mean 
temperature  value  of  the  longitudinal  acoustic  as  well  as  of  optical  phonons  taking  place  in  the 
interaction  [2-3,  10-12].  The  average  phonon  temperature  0  is  related  to  the  average  phonon 
energy  EP  via  the  Boltzmann  constant  kB:  ©=  EP  /  kB.  The  model  presented  in  Equation  (1) 
exhibits  a  good  fit  to  our  PL  experimental  data  of  the  GaN  film,  as  can  be  seen  in  Figure  2.  The 
fit  resulted  in  the  following  values:  E(0)=3.4866  eV,  a=55  meV,  and  0=309  K  and  thus  EP  =27 
meV.  These  values  agree  with  the  ones  obtained  previously  for  GaN  films  via  the  same  model 
utilizing  optical  absorption  [13]  and  contactless  electroreflectance  [10]. 

As  can  be  observed  in  Figure  2  the  temperature  dependence  of  the  PL  emission  of  the 
GaN/Alo^Gao.gN  superlattice  does  not  follow  the  monotonic  variation  of  Equation  1  for  the  band 
edge  to  lower  energies  with  increasing  temperature.  Similar  behavior  has  previously  been 
observed  in  the  quantum  well  structures  of  AlGaN/InGaN/GaN  [4],  InGaAs-InP  [5],  and 
GalnAs/AlInAs  [6]  as  well  as  in  superlattices  of  ZnSe/MgS  [7]  and  GaAs/AlAs  [8].  This 
anomalous  emission  characteristic  has  been  attributed  to  the  presence  of  density  of  states  (DOS) 
which  is  due  to  a  certain  degree  of  disorder  occurring  mainly  at  interfaces,  which  may  be  of 
compositional  and/or  structural  origin  [5].  An  explanation  of  the  anomalous  behavior  of  the 
excitonic  luminescence  was  given  in  [5]  in  terms  of  this  disorder;  we  apply  this  explanation  to 
our  results  presented  in  Figure  2.  As  the  temperature  rises  from  10K  to  50  K,  the  excitons  gain 
just  enough  thermal  energy  to  become  trapped  at  the  lower  energy  levels  of  the  DOS,  which 
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results  in  the  red  shift  of  the  PL  emission.  From  50K,  a  thermal  equilibrium  is  achieved  and  the 
excitons  populate  the  higher  levels  of  the  DOS  (blue  shift)  up  to  a  maximum  energy  level,  in  this 
case  ~  130K  corresponding  to  an  emission  energy  of  3.8  eV.  From  that  point  on,  a  “normal” 
band  gap  emission  takes  over  which  can  be  described  in  terms  of  Equation  1. 

A  model  has  been  recently  proposed  to  explain  the  blue  temperature-induced  shift  in 
InGaN-based  quantum  well  luminescence  [4],  The  model  is  based  on  band-tail-  filling  of  a 
Gaussian  DOS  with  a  parameter  a  which  describes  the  dispersion  of  the  DOS  (i.e.  its  width). 
The  model  calculation  is  given  by  [4] 

2 

E(T)=ED-J-  (2) 

where  Ed  is  attributed  to  the  energy  range  between  the  center  of  the  Gaussian  DOS  of  the 
electrons  and  that  of  the  holes.  In  addition,  the  model  assumes  that  ED  depends  on  temperature  in 
the  same  manner  as  the  gap-energy,  and  thus  ED  can  be  described  via  Equation  1.  This  model  is 
valid  in  the  temperature  range  where  thermal  equilibrium  has  been  established:  in  our  case  -  50K 
and  above.  We  use  this  model  to  analyze  our  data  in  order  to  estimate  the  extent  of  the  DOS  in 
the  superlattice;  the  fit  is  presented  in  Figure  2.  We  obtain  the  value  for  a  to  be  ~  11  meV, 
which  is  a  relatively  small  compared  to  the  value  of  the  superlattice  emission  (~  3.8  eV) 
indicating  the  high  quality  of  the  SL  interfaces.  Additionally,  the  inset  to  Figure  2  shows  the  PL 
linewidth  as  a  function  of  temperature.  The  linewidth  characteristic  can  also  be  separated  into 
three  different  temperature  regimes  similar  to  those  observed  for  the  PL  emission;  the 
broadening  mechanisms  will  be  addressed  at  a  later  study.  Equation  2  yields  the  values  of  the 
average  electron-phonon  interaction  strength  a=  135  meV  and  of  the  mean  phonon  temperature 
0  =  573  K,  the  latter  which  corresponds  to  a  mean  phonon  energy  Ep=50  meV.  This  higher  value 
of  Ep  obtained  for  the  SL  relative  to  that  of  GaN  (Ep=27meV)  may  be  explained  in  view  of  the 
model  calculations  developed  for  polar-optical  phonons  in  GaN/AlGaN-based  quantum  wells  of 
wurtzite  structure  [14-15].  In  that  study,  Komirenko  et  al.  have  shown  that  phonon  modes  in  the 
heterolayers  of  GaN-AlGaN  based  wurtzite  heterostructures  can  exhibit  the  higher  energies  of 
the  modes  of  the  AlGaN  material  [14,  15],  In  light  of  the  above  model,  the  higher  average 
phonon  energy,  Ep,  found  in  our  GaN  wells  are  expected  based  on  the  predicted  phonon  mode 
structure  in  GaNM.l0.2Gao.gN  barriers. 
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In  order  to  gain  further  insight  into  the  mechanisms  determining  the  recombination 
dynamics  of  the  superlattice,  the  luminescence  intensity,  I(T),  as  a  function  of  temperature  was 
investigated.  Figure  3  presents  the  experimental  data  of  that  study  and  the  curve  fit  to  the  data 
obtained  from  the  relation  which  describes  the  PL  intensity  in  the  presence  of  a  non-radiative 
channel  with  thermal  activation  energy  EA[16]: 

Ml _ L _  0) 

I0  l  +  Cexp(-EA/kBT ) 

In  this  relation  C  is  a  constant  and  I0  is  the  PL  intensity  at  zero  absolute  temperature.  An 
activation  energy  EA=16  meV  was  ascertained  from  the  fit.  This  value  is  in  the  same  order  of 
magnitude  and  agrees  with  that  obtained  for  the  energy  width  of  the  DOS  (11  meV),  the  small 
difference  being  attributed  to  an  experimental  error  and  the  approximate  nature  of  the  above 
models.  Thus  at  higher  temperatures  the  tail  states,  due  to  interface  disorder,  are  the  non- 
radiative  centers  in  the  superlattice  emission. 

The  temperature  dependence  of  the  luminescence  blue  shift  in  InGaN  based  quantum 
wells  has  been  discussed  in  terms  of  the  screening  of  the  piezoelectric  field  [17].  In  that  model, 
as  the  temperature  increases  free  carriers  are  thermally  activated  which  in  turn  causes  the 
screening  of  the  field.  As  a  result,  modification  of  the  band-gap  takes  place  that  manifests  itself 
in  the  blue  shift  of  the  emission.  Although  this  effect  may  be  a  possible  additional  mechanism 
causing  the  temperature  dependent  blue  shift  of  the  luminescence,  it  does  not  explain  the  initial 
red  shift  of  the  emission  energy.  As  mentioned  previously,  similar  luminescence  behavior  has 
been  observed  in  the  zincblende  IH-V  quantum  wells  [5-8]  which  have  considerably  smaller 
piezoelectric  fields  [18].  We  do  not  however  preclude  the  screening  effect  of  the  PEF  as  an 
additional  mechanism  for  the  blue  shift  of  the  PL 

In  conclusion,  the  temperature  behavior  of  the  luminescence  in  a  GaN/Alo^GaosN 
superlattice  was  investigated.  Our  analysis  verifies  the  existence  of  a  defect  density  of  states,  in 
the  recombination  of  excitons  at  low  temperatures.  We  also  find  that  at  high  temperatures,  the 
density  of  states  acts  as  a  non-radiative  channel  for  the  luminescence.  The  higher  value  of  the 
average  energy  of  the  phonons  involved  in  the  superlattice  scattering  is  consistent  with  the 
predictions  for  the  phonon  modes  of  GaN/AlGaN  heterolayers. 
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FIGURE  CAPTIONS 


1.  The  PL  at  lOKof  the  SL.  The  peak  at  3.795  is  due  the  emission  from  the  GaN  wells.  The 
inset  to  the  Figure  depicts  the  PL  intensity  as  a  function  of  the  laser  excitation  power. 

2.  The  temperature  dependence  of  the  PL  peak  position  of  the  SL  GaN  wells  (circles)  and  of 
GaN  film  (squares).  The  line  represent  the  model  of  Equation  2.  The  inset  to  the  Figure 
depicts  the  PL  linewidth  as  a  function  of  temperature. 

3.  The  normalized  PL  intensity  of  the  SL  GaN  wells  as  a  function  of  temperature;  Data  points 
and  the  curve  fit  (line)  to  Equation  3. 
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ABSTRACT 

This  paper  presents  a  study  of  the  quasi-LO  and  -TO  modes  in  Wurtzite  AIN  which  originate  from 
the  interaction  of  phonons  belonging  to  the  A1  and  El  symmetry  groups.  In  order  to  analyze  the 
allowed  quasi  as  well  as  pure  Raman  modes,  the  modes  were  observed  in  a  rotating 
crystallographic  coordinate  system,  and  the  Raman  tensors  of  the  Wurtzite  crystal  structure  were 
calculated  as  a  function  of  the  crystallographic  rotation.  The  frequencies  of  the  quasi-modes  of 
Wurtzite  AIN  were  also  analyzed  in  terms  of  the  interaction  of  the  polar-phonons  with  the  long 
range  electrostatic  field  model.  The  experimental  values  of  the  Raman  frequencies  of  the  quasi¬ 
phonons  concur  with  these  expected  from  the  model,  implying  that  the  long  range  electrostatic  field 
dominates  the  short  range  forces  for  polar  phonons  in  AIN. 


INTRODUCTION 

Raman  spectroscopy  has  proven  to  be  a  useful  and  informative  tool  in  the  investigation  of 
semiconductors  with  the  Wurtzite  (WZ)  crystal  structure,  in  particular  in  the  investigation  of 
material  quality  as  well  as  in  the  study  of  phonon  dynamics  and  interactions  [1-9].  Due  to  the  long 
range  electrostatic  field  associated  with  the  LO  phonons,  various  types  of  interactions  are  allowed 
in  the  WZ  semiconductors,  among  them  the  phonon-plasmon  coupling  [3].  Another  is  the 
possibility  of  interaction  between  polar  modes  belonging  to  different  symmetry  groups  which 
result  in  the  creation  of  new  mixed-symmetry  phonons  referred  to  as  quasi  modes  [1],  This  paper 
presents  a  study  of  the  LO  and  TO  quasi  modes  in  AIN  crystallite  which  originate  from  the 
interaction  of  phonons  belonging  to  the  A1  and  El  symmetry  groups;  the  study  focuses  on  the 
mode-selection  rules  and  mode-frequency  calculations  as  well  as  on  the  Raman  technique  which 
enables  the  observation  of  the  modes. 

The  theory  developed  by  Loudon  formulates  that  in  uniaxial  materials  the  polar  phonon 
characteristics  may  be  affected  via  two  interaction  mechanisms:  one  due  to  the  long  range 
electrostatic  field,  and  the  other  due  to  the  short  range  field  which  exhibits  the  anisotropy  of  the 
vibrational  force  constants  [10].  The  phonon  dynamics  and  thus  the  Raman  spectra  depend  on 
which  of  the  two  mechanism  is  the  dominant  interaction. 

For  the  case  where  the  long  range  electrostatic  field  is  the  dominant  mechanism,  the 
interaction  of  the  polar  phonons  with  the  long  range  electrostatic  field  may  result  in  a  significant 
frequency  separation  between  the  group  of  the  TO  phonons  relative  to  that  of  the  LO  phonons,  as 
is  schematically  depicted  in  Figure  la.  Moreover,  the  TO  phonons  belonging  to  different 
symmetry  are  grouped  together  in  a  relatively  narrow  frequency  range;  the  same  holds  for  the  LO 
phonons.  One  consequence  of  the  dominant  electrostatic  field  interaction  is  that  the  phonons 
exhibit  nearly  pure  longitudinal  or  transverse  character:  essentially  no  mixing  occurs  between  the 
transverse  and  the  longitudinal  atomic  vibrations.  The  other  consequence  is  that  under  certain 
propagation  and  polarization  conditions,  phonons  of  mixed  A1  and  El  symmetry  character  exist 
and  can  be  observed  in  the  Raman  spectra.  These  mixed  symmetry  modes  are  termed  quasi-LO 
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and  -TO  modes.  The  frequencies  of  such  quasi  modes  are  predicted  by  the  theory  to  be  between 
the  values  of  the  pure-Al  and  the  pure-El  for  each  of  the  LO  and  the  TO  bands. 

Alternatively,  for  the  case  where  the  short  range  interatomic  forces  are  dominant,  the  LO- 
TO  splitting  will  be  small  and  in  this  case  the  TO  and  LO  mode  of  each  symmetry  group  will  occur 
in  a  relatively  narrow  band  as  shown  in  Figure  lb.  The  theory  predicts  that  the  phonons  will 
exhibit  essentially  pure  symmetry  character  (either  A 1  or  El).  A  consequence  is  that  under  certain 
propagation  conditions  the  phonons  will  exhibit  mixing  of  the  transverse  and  longitudinal 
vibrations  for  a  given  symmetry.  These  mixed  LO-TO  modes  are  termed  quasi-El  and  -A1  modes. 
The  frequencies  of  such  quasi  modes  are  predicted  by  the  theory  to  be  in  bands  between  the  values 
of  the  TO-LO  modes  of  the  pure  A1  or  El  phonons. 

In  AIN  the  El(TO)  and  Al(TO)  Raman  frequencies  are  grouped  together  in  ~  60  cm'1 

frequency  range,  the  El(LO)  Al(LO)  in  ~  20  cm’1  range,  and  the  LO-TO  group-splitting  is  ~  220 

cm'1.  The  frequency  scheme  of  AIN  thus  implies  the  dominance  of  the  long  range  electrostatic 

force  interaction.  The  mode-mixing  in  AIN  may  occur  if  the  propagation  direction  (q-vector)  of  the 
quasi-polar  phonons,  their  polarization  LO,  and  one  of  their  two  TO’s,  lie  in  the  plane  which  spans 
between  the  c  and  the  a,  (or  aj)  crystallographic  axes  [1].  This  plane  is  referred  as  to  the  mixing 
plane.  For  example,  the  pure-Al  phonon  has  a  c-direction  polarization  while  the  pure-El  phonon 
is  polarized  in  the  basal  plane;  thus  a  quasi-phonon  with  the  q-vector  between  the  c  and  the  a,  axes 
would  exhibit  a  mixed  polarization  of  Al-El  symmetry.  When  the  q-vector  lies  along  the 

crystallographic  axes  or  in  the  ara2  plane  only  pure  phonons  are  observed  in  the  spectra. 

Arguello  et  al.  investigated  the  long  range  electrostatic  field  interaction  and  the  quasi  modes 
in  the  WZ  materials:  ZnO,  ZnS,  and  CdS  [1].  In  their  study  a  frequency  shift  was  observed  of  the 
quasi  modes  relative  to  that  of  the  pure  modes,  in  agreement  with  Loudons’  theory.  Filippidis  et. 
al.  have  carried  out  plane-wave  pseudopotential  calculations  to  investigate  the  optical  modes  in 
GaN  and  AIN;  their  model  and  experiment  concerning  the  Al(TO)  mode  agree  with  that  of 
Loudons’  [9]. 
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In  our  study  we  investigate  the  behavior  of  the  quasi-TO  and  the  quasi-LO  in  AIN 
crystallite.  We  present  calculations  of  the  WZ  Raman  tensors  as  a  function  of  rotation  angle  from 
which  the  allowed  modes,  pure  and  quasi,  can  be  deduced  and  analyzed  for  different  crystal 
orientations.  The  experimental  values  of  the  quasi  frequencies  agree  well  with  those  predicted  by 
the  phonon-electrostatic  field  interaction  model.  The  motivation  behind  the  present  study  was  to 
gain  information  on  the  Raman  frequencies  from  scattering  geometries  which  are  not  along  the 
crystallographic  directions.  As  mentioned  above,  Raman  spectroscopy  has  been  successfully 
utilized  in  the  characterization  of  the  group  Ill-nitrides  semiconductors  and  has  proven  useful  in 
studies  of  the  stress  state  of  films  grown  on  the  c-axis.  However  for  films  grown  on  other 
crystallographic  directions,  a  technique  which  potentially  might  minimize  film  substrate  interfacial 
stress,  the  contribution  of  the  long  range  electrostatic  field  to  the  Raman  shift  has  to  be  taken  into 
consideration.  This  contribution  has  to  be  deconvoluted  it  in  order  to  isolate  the  Raman  shift  due  to 
stress. 

EXPERIMENT 

The  micro-Raman  scattering  experiments  were  carried  out  at  room  temperature  utilizing  the 
514.5  nm  line  of  the  ion  argon  laser  and  the  J-Y  U1000  scanning  double  monochromator.  The  spot 

size  of  the  laser  on  the  sample  was  -2  pm  in  diameter.  The  experimental  error  of  the  Raman  peak 

position  and  the  angle  of  rotation,  0,  is  ±0.5  cm'1,  and  ±5°  respectively.  The  AIN  crystallite, 

depicted  in  Figure  2,  was  grown  via  the  sublimation  process  and  its  size  is  -500x25x25  |i.m  [5]. 

In  the  figure,  the  face  of  the  (lOlo)  plane  is  along  the  c-axis  while  the  face  of  the  (0001)  plane  is  at 
the  base  of  the  hexagonal  crystallite  [5]. 
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RESULTS  AND  DISCUSSION 

In  this  study  we  observed  the  AIN  quasi-TO  and  -LO  modes  by  rotating  the  crystallite 
relative  to  the  laboratory  coordinate  system  as  schematically  depicted  in  Figure  3.  In  the  figure, 
the  X,  Y,  and  Z  are  the  axes  of  the  fixed  coordinate  laboratory  system  (where  X  is  out  of  the  paper 
plane),  while  the  c,  a,  and  %  are  the  axes  of  the  rotating  coordinate  system  of  the  crystal  (where  a2 
is  out  of  the  paper  plane  and  is  coincided  with  the  X-axis),  and  N  is  the  normal  to  the  crystallite 

surface  (and  is  coincided  with  the  a,-axis).  The  AIN  crystallite  was  rotated  through  an  angle  0 

about  the  X-axis  such  that  the  quasi-phonon  propagation  direction,  the  q-vector  (q),  as  well  as  its 
LO  and  one  of  its  two  TO  components  lay  in  the  mixing  c-a,  plane  (the  other  TO  is  along  the  X- 
axis).  The  micro-Raman  spectra  was  acquired  in  a  back  scattering  geometry  from  the  Y-axis  for 

various  values  of  0. 

In  order  to  investigate  the  polarization  conditions  upon  which  the  quasi  modes  appear  in  the 
spectra,  the  Raman  tensors,  Re,  of  the  rotating  crystal  (relative  to  a  fixed  X,Y,Z  laboratory 

coordinate  system)  are  calculated  via  the  rotational  transformation  equation  [11]: 

Re  =  txrtx ■'  (i) 


where  Tx  is  the  rotation  matrix  about  the  X-axis  [12].  The  Raman  tensors  of  the  WZ  structure  take 
the  form: 
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The  tensors  describe  the  in-crystal  projections,  due  to  the  rotation,  of  the  incident  and  scattered 
polarizations  of  the  light.  The  group  theory  analysis  indicates  that  the  A1  is  a  non  degenerate 
mode  and  thus  has  a  one  matrix  representation,  while  the  El  and  E2  are  both  doubly  degenerate 
modes  and  therefore  each  mode-representation  requires  two  matrices.  The  entries  of  each  matrix, 
i.e.  the  constants  h,  b,  d,  and  f,  are  related  to  the  Raman  cross  section  which  determines  the 
Raman  intensity  of  each  of  the  modes  [13].  Moreover,  for  each  polar  mode  the  direction  of  the 
polarization  is  indicated  in  the  subscripts:  A1  is  polarized  along  the  c  axis  while  El  is  polarized 
along  a,  and  a,,  the  latter  which  coincides  with  the  X-axis  (Elx).  The  polarization  direction  of  the 
polar  modes,  like  the  Raman  tensors,  can  be  expressed  in  terms  of  the  laboratory  coordinate 

system;  however,  it  is  not  necessary  for  the  present  study.  For  0=0  the  crystal  coordinate  system 

coincides  with  that  of  the  laboratory  (a,=Y  and  c=Z,  see  Fig.  3),  and  Eq.  2  takes  the  usual  form  [1, 
13]. 

For  the  first  set  of  experiments  we  chose  the  Raman  setup  such  that  the  incident  light  is 
polarized  along  the  X-axis  and  the  scattered  light  has  both  X  and  Z  polarization  directions, 

expressed  in  Raman  notation  as:  Y^X  ^Y.  This  ch°*ce  of  polarization  implies  that  the  XX  and 

the  XZ  components  of  each  of  the  Raman  tensors  will  contribute  a  Raman  signal  provided  they  are 
not  zero.  Figure  4  presents  the  Raman  spectra  for  this  setup  for  rotations:  0  =  0°,  35°,  and  70°. 
Inspection  of  Eq.2  indicates  that  the  E2  mode  should  appear  in  the  spectra  since  the  XX  and  the 
XZ  components  of  its  tensors  are  not  zero;  our  spectra  display  the  E2  mode  at  655  cm”1.  This 

mode  is  non-polar  and  as  such  no  change  in  frequency  is  expected  as  a  function  of  the  rotation. 
The  XZ  component  of  the  Elx  is  non-zero,  and  although  the  q-vector  of  the  mode  lies  in  the  a-c 
mixing  plane  its  polarization  is  not:  it  is  in  the  X-direction  (TO  component)  and  thus  no  symmetry 
mixing  is  expected.  This  analysis  indicates  that  the  spectra  should  exhibit  the  pure-El(TO)  with  no 
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change  of  frequency  upon  rotation.  The  pure-El(TO)  in  our  spectra  is  at  669  cm'1  and  is 
independent  of  rotation. 

Next  we  investigated  the  A1  mode  which  contributes  a  signal  to  the  Raman  spectra  via  the 
non-zero  value  of  its  XX  tensor  component.  The  pure-Al  mode  has  solely  c-polarization; 
however,  upon  rotating  the  crystal  quasi  modes  are  created  for  which  their  TO  and  LO  components 
have  a  c-polarization  component  of  A1  symmetry  as  well  as  an  a, -polarization  component  of  El 
symmetry.  This  situation  is  schematically  depicted  in  the  inset  to  Figure  3.  In  the  figure,  the  TO 
component  of  the  quasi-mode  lies  in  the  mixing  plane  and  is  composed  of  A1  and  El  polarization 
symmetry  components;  the  quasi-TO  mode  will  have  mostly  A1  symmetry  for  q-vectors  near  the 
a-axis  and  mostly  El  symmetry  for  q-vectors  near  the  c-axis.  The  frequency  of  the  quasi-TO, 
according  to  Loudons’  theory,  should  be  between  the  frequency  of  the  pure-Al  (TO)  mode  and  that 
of  the  pure-El(TO)  mode.  As  can  be  seen  in  Figure  4,  the  frequency  of  the  Al(TO)  mode  shifts 
from  its  pure  value  at  608  cm'1  towards  that  of  the  El  (TO)  mode  as  a  function  of  the  rotation  angle: 

for  0=35°  the  peak  of  the  quasi-TO  is  at  614  cm'1  and  for  0=70°  is  at  620  cm'1.  The  spectrum  of 

0=70°  exhibits  a  small  peak  of  the  pure-Al(TO)  and  is  attributed  to  an  internal  reflection. 

The  frequencies  of  the  quasi-TO  can  be  calculated  via  the  relation  [1,10] 

^Q(TO)  =  ®El(TO)COS  (90-P)  +  <^Al(TO)S^n  (90  -  [3)  (3) 

where  p  (see  Fig.  3)  is  the  propagation  angle  of  the  q-vector  of  the  quasi-phonon  (q)  in  the 

crystallite,  and  is  related  to  the  angle  of  rotation,  0,  via  Snell’s  law  and  the  momentum- 

conservation  law  of  the  Raman  scattering  effect  as  is  discussed  next.  The  conservation  law  states 
that  in  the  scattering  medium  q  =  k,  -  ks  where  k,  and  ks  are  wavevectors  of  the  incident  and 
scattered  photons  respectively  [13-14],  Thus,  for  our  backscattering  geometry  the  q-vector  lies 

along  the  same  path  as  that  of  the  incident  and  scattered  photons.  To  find  the  angle,  P,  of  that  path 
the  index  of  refraction  has  to  be  taken  into  account  via  Snell’s  law:  n,sin0=n2sinP,  where  for  air 
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n,=l  and  for  AIN  n2=2.2  which  is  the  accepted  value  in  the  visible  range  and  of  negligible 
anisotropy  [15].  Our  calculations  yield  a  value  of  p=15°  for  the  rotation  0=35°,  and  from  Eq.  3, 

ooQ(TO)  =  612  cm'1,  a  value  which  concurs  with  our  experimental  finding  of  614  cm’1.  For  the  7Cf 

rotation,  (3=25°  and  C0q(TO)  =  619  cm'1,  a  value  which  is  similar  to  the  experimental  result  of  620 

-l 

cm  . 

Lastly  we  investigated  the  quasi-LO  Raman  mode.  In  order  to  observe  this  mode  the 
incoming  polarization  of  the  light  was  set  at  the  Z-direction  such  that  the  Raman  geometry  was 

Y(z  ^ .  This  Raman  geometry  enables  the  observation  of  the  quasi-LO  modes  since  the  ZZ 

component  of  the  Ela  Raman  tensor  is  non-zero  for  9>0  as  is  indicated  by  Eq.  2.  In  this 

configuration  the  q-vector  of  the  phonon  as  well  as  its  polarization  directions  lie  in  the  mixing 
plane.  Figure  5  presents  the  Raman  spectra  of  the  quasi-LO  mode  for  0=15°,  35°,  and  70°,  and  the 

pure  Al(LO)  mode  at  0=90°;  the  Raman  peak  positions  are  at  913  cm'1,  912  cm'1,  908  cm'1,  and 

890  cm'1  respectively.  Taking  into  account  the  index  of  refraction  and  utilizing  the  equation  for  the 
quasi-LO  frequencies  [1,10], 

^QtLO)  =  ®Al(LO)COS  (^“P)  +  tt)El(LO)S^n  (^0-P) 

where  C0A1(LO)  is  890  cm'1  and  (t>E1(Lc»  is  916  cm'1  [4],  we  found  the  calculated  quasi  frequencies  to 

be  915  cm'1,  914  cm'1,  and  91 1  cm'1,  which  agree  well  with  the  above  experimental  values. 

Figure  6  summarizes  our  findings.  In  the  figure  the  lines  denote  the  TO  and  LO  quasi¬ 
frequencies  calculated  utilizing  Loudons’  model  (Eqs.  3  and  4),  and  the  dots  represent  our 
experimental  values  for  the  modes.  The  additional  data  points  of  the  quasi-TO  modes,  represented 
by  X’s,  are  from  Filippidis  et  al  who  were  the  first  to  observe  the  quasi-TO  modes  in  AIN  [9];  our 
results  are  in  agreement  with  theirs.  Moreover,  as  depicted  in  the  figure,  our  study  of  the  quasi-TO 
as  well  as  the  quasi-LO  concur  with  Loudons’  model,  implying  the  presence  of  a  long  range 
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electrostatic  field  which  is  active  in  the  coupling  of  modes  belonging  to  the  different  symmetry 
groups. 

One  important  implication  of  the  above  findings  involves  the  accurate  determination  of  the 
Raman  frequency  acquired  from  AIN  thin  films.  In  such  films,  due  to  their  micron-size  dimension 
along  one  of  the  crystallographic  axes,  the  determination  of  the  various  Raman  frequencies  requires 
the  data  to  be  acquired  at  grazing  angle  geometry.  Such  an  experimental  setup,  which  deviates  by 
an  angle  from  the  crystallographic  coordinate  system,  may  result  in  a  mode  mixing  of  a 
characteristic  Raman  frequency.  Due  to  the  strong  angular  dispersion  of  the  quasi-TO  mode 
frequency,  a  disparity  in  the  TO  Raman  frequency  values  of  thin  films  is  expected  to  be  observed; 
less  disparity  is  expected  for  the  LO  mode.  A  similar  effect  in  which  the  frequency  is  shifted  from 
its  pure  value  is  also  expected  to  occur  in  films  grown  on  crystallographic  directions  which  do  not 
coincide  with  the  c-axis.  As  Raman  spectroscopy  is  widely  used  in  the  determination  of  stress  in 
thin  films  via  the  Raman  frequency  shift,  the  possible  additional  shift  due  to  the  electrostatic  field 
effect  has  to  be  considered. 

CONCLUSIONS 

The  Raman  tensors  of  the  WZ  structure  were  calculated  as  a  function  of  the  crystallographic 
rotation  and  the  allowed  modes  for  various  rotations  were  studied.  The  LO  and  TO  quasi-modes 
of  AIN  crystallite  were  observed  and  analyzed  in  terms  of  the  phonon-electrostatic  field  interaction 
model.  The  Raman  selection  rules  as  well  as  the  experimental  values  of  the  Raman  frequencies 
concur  with  the  model  implying  the  dominance  of  the  long  range  electrostatic  field  in  AIN.  The 

quasi-TO  mode  exhibits  a  shift  of  12  cm'1  relative  to  its  pure  value  for  mode-propagation  angle  ~ 
25°  measured  from  the  a-axis.  For  the  same  propagation  direction,  the  quasi-LO  mode  exhibits  a 
shift  of  -  5  cm'1. 
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FIGURE  CAPTIONS 


Figure  1.  Raman  frequency  scheme  and  directions  of  phonon  vibration  for  (a)  the  case  when  the 
long  range  electrostatic  field  is  the  dominant  interaction,  and  (b)  when  the  anisotropy  of  the  short 
range  dominates  the  interaction. 


Figure  2.  A  micrograph  of  the  Wurtzite  AIN  crystallite,  size  -  500x25x25  pm.  Along  the  c-axis 

is  the  face  of  the  (lOlo)  plane,  the  face  of  the  (0001)  plane  is  at  the  base  of  the  hexagonal 
crystallite. 


Figure  3.  The  geometry  of  the  experimental  setup.  The  laboratory  coordinate  system  is:  X,  Y, 
and  Z  (X  is  out  of  the  paper  plane).  The  crystal-rotating  system  is:  c,  a,,  and  a-,  (%  out  of  the  paper 
plane  and  is  coincided  with  the  X-axis).  N  is  the  normal  to  the  surface  and  coincides  with  the  a,- 

axis.  The  angle  of  rotation  about  the  X-axis  is  0,  while  the  angle  of  refraction  and  phonon 

propagation  is  [3 .  The  dashed  line  represents  the  propagation  of  the  phonon  with  a  q- vector;  the 

arrows  represent  its  LO  and  TO  components.  The  Raman  spectra  were  acquired  in  a  backscattering 
geometry  from  the  Y.  The  inset  to  the  figure  describes  the  Al-El  mixing  of  the  LO  and  the  TO 
components. 


Figure  4.  The  Raman  spectra  for  rotation  angle  0=  0°,  35°,  and  70°.  The  pure-Al(TO)  is  at  608 

cm'1,  and  the  quasi-TO  modes  are  at  614  cm'1,  and  620  cm'1.  All  intensities  are  normalized  to  that 
of  the  E2  mode. 
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Figure  5.  The  normalized  Raman  spectra  of  the  quasi-LO  modes  for  rotation  angle  0=  15°,  35°, 


and  70°.  The  pure  Al(LO)  is  also  presented  for  which  0  =90°. 


Figure  6.  The  behavior  of  the  quasi-LO  and  -TO  modes  as  a  function  of  phonon  propagation 
angle  (3  (measured  from  the  a-axis).  The  lines  represent  the  TO  and  LO  quasi  frequencies 

calculated  using  Loudons’  model  equations  3  and  4;  the  dots  represent  our  experimental  values  for 
the  LO  and  the  TO  modes,  the  value  of  the  pure  El(LO)  is  from  reference  4  (presented  by  a 
triangle),  while  the  X’s  represent  the  data  points  of  the  quasi  TO  from  reference  9. 
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Abstract  This  paper  reviews  Raman  studies  of  GaN,  AIN,  InN,  and  the  alloy  systems  Ga,Al,J\T 
and  In.xGai.jtN-  The  review  focuses  on  the  applications  of  Raman  spectroscopy  to  material 
characterizations  and  phonon  dynamics  of  the  wurtzite  (WZ)  and  the  zincblende  (ZB)  polytypes. 
Among  the  topics  addressed  are  structure  determination  and  microstructure  identification,  stress 
and  strain  in  film-substrate  interfaces  and  in  crystallites,  phonon-plasmon  coupling,  and  the 
determination  of  free  carrier  density.  The  issue  of  phonon  mode-type  in  the  alloys  is  also 
examined,  as  well  as  the  topic  of  alloy  order.  The  subject  of  the  quasi-modes,  their  angular 
dispersion,  and  the  dominance  of  the  long  range  electrostatic  field  is  addressed.  Also,  the  issues  of 
the  isotope  effects  on  the  various  vibrational  modes  as  well  as  the  phonon  lifetimes  are  discussed. 
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1.  Introduction 


The  potential  of  UV  optoelectronic  devices  made  of  III-V  nitrides  has  prompted  numerous  studies 

i 

focusing  on  the  properties  of  this  family  of  materials  [1-5].  One  of  the  most  attractive  properties 
of  the  III-V  nitrides  is  the  nature  of  their  wide  band  gaps,  which  in  the  wurtzite  structure  have  the 
values  of  1.95  eV  for  InN,  3.45  eV  for  GaN,  and  6.4  eV  for  AIN  [2,  5].  In  addition,  the  superior 
mechanical  and  thermal  properties  of  the  III-V  nitrides,  as  well  as  the  feasibility  of  alloying  and 
the  material  polymorphism,  have  made  this  family  of  materials  a  unique  candidate  for  synthesizing 
high-performance  UV  devices. 

Raman  spectroscopy  has  contributed  a  great  deal  to  the  advances  in  the  III-V  nitride  field,  and 
the  key  studies  leading  to  these  advances  are  reviewed  here.  This  review  focuses  in  particular  on 
the  application  of  Raman  spectroscopy  to  topics  that  are  of  concern  to  material  scientists 
including  structure  determination,  stress  analysis,  determination  of  free  earner  concentrations,  and 
material  quality  characterization.  Other  issues  such  as  alloy-disorder  and  alloy-mode,  isotopic 
effects,  and  aspects  of  phonon  dynamics  and  phonon  lifetimes  are  also  addressed.  The 
implications  of  these  issues  for  material  characterizations  are  discussed  as  well. 

In  Section  2  a  discussion  of  the  underlying  principles  involved  in  Raman  scattering 
experiments  on  wide  band-gap  semiconductors  is  presented.  Section  3  focuses  on  Raman 
selection  rules  and  studies  concerning  the  structural  properties  of  the  III-V  nitrides.  In  particular, 
Section  3.1  surveys  polarized  Raman  analyses  of  the  wurtzite  and  zincblende  structures,  and 
Sections  3.2  and  3.3  review  structure-related  studies  of  GaN,  and  AIN  and  InN  respectively. 
Stress  effect  studies  on  GaN  and  AIN  films  and  crystallites  are  next  reviewed  in  Section  4.  The 
relevance  of  the  hydrostatic  and  biaxial  Raman  stress  coefficients  to  stress  analysis  of  the  nitride 
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films  and  crystals  is  discussed.  Specifically,  Sections  4.1  and  4.2  focus  on  stress  analysis  of  GaN 
and  AIN,  respectively. 

Section  5  is  devoted  to  Raman  analyses  of  the  LO-  and  7'0-quasi-modes  and  their  frequency 
dispersion;  the  implications  of  the  dispersion  to  film  characterization  are  also  presented.  Another 
topic  addressed  in  Section  6  is  the  LO  phonon-plasmon  coupling  and  the  resulting  spectral 
lineshape.  The  validity  of  using  the  lineshape  as  a  measure  of  the  free  carrier  concentration  in  GaN 
is  discussed.  In  Section  7  studies  of  the  isotopic  effects  on  the  phonon  dynamics  as  well  as  studies 
concerning  the  mode  lifetimes  are  summarized.  Section  8  presents  Raman  studies  of  III-V  nitride 
alloy  systems.  The  one-mode  and  two-mode  behavior  of  phonons  in  the  alloy  is  described,  and 
experimental  and  theoretical  studies  are  presented.  The  topic  of  alloy  disorder  and  its  detectability 
via  Raman  spectroscopy  is  addressed  as  well.  Section  9  concludes  with  comments  about 
directions  of  future  research. 


2.  Experimental  Considerations  for  Raman  Scattering  of  Wide  Band-Gap  Semiconductors 

In  its  general  form  the  Raman  scattering  intensity  1  can  be  expressed  as 


I(coL)x  6f||es  •  R 


-  |2 


%(Ea  -h(DL){Ep-tUDs) 


(1) 


where  ox,  ah  are  the  incoming  and  scattered  laser  frequencies,  respectively;  Ea  and  Ep  are  the 
energies  of  intermediate  crystal  states  (to  be  defined);  R  is  the  Raman  tensor;  and  esand  e L  are 
the  scattered  and  incident  polarization  vectors.  In  (1)  the  first  term  is  due  to  the  dipole  transition 
radiation,  the  second  represents  the  Raman  selection  rules,  which  come  about  from  crystal 
symmetry  considerations,  and  the  last  term  leads  to  resonance  effects  [6]. 


Energy  conservation  relates  ox  and  ah  to  the  phonon  frequency  ay. 
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As  the  incident  laser  frequency  approaches  the  energies  of  the  crystal  intermediate  states,  the 
Raman  intensity  becomes  larger  and  the  signal  is  said  to  be  resonance-enhanced,  an  additional 
enhancement  comes  from  the  dipole  radiation  term.  In  general,  for  semiconductors  there  may  be 
three  relevant  types  of  intermediate  states:  Bloch  states  which  are  the  conduction-valence  bands, 
exciton  states,  and  in-gap  impurity  states.  Thus  as  the  incident  frequency  approaches  the  band-gap 
frequency,  an  enhancement  of  the  Raman  intensity  should  be  observed.  Similarly,  resonance  may 
be  achieved  via  the  interaction  of  the  incoming  light  with  the  exciton  states  as  well  as  the  impurity 
states. 

The  band  gaps  of  InN,  GaN,  and  AIN  are  listed  in  Table  1  [2].  It  is  evident  that  if 
enhancement  of  the  Raman  signal  of  pure  GaN  is  to  be  achieved,  a  laser  line  in  the  lower  energy 
range  of  the  UV  (*  2.7—4  eV)  should  be  utilized  [7],  unless  impurity  resonance  is  present  [8]. 
Resonance  Raman  scattering  of  AIN  and  AlGaN  alloys,  on  the  other  hand,  requires  laser 
excitation  lines  in  the  deep  UV  range  (*  4-7  eV)  depending  on  the  alloy  composition  [9].  The 
resonance  of  the  Raman  signal  of  InN,  due  to  its  narrower  band  gap,  can  be  achieved  via 
excitation  in  the  visible  optical  range.  For  laser  excitation  energy  much  above  the  fundamental 
band-gap  energy  of  the  semiconductor,  the  light  is  absorbed  and  thus  only  a  small  volume  of  the 
material  is  probed.  In  the  case  of  strong  absorption  the  Raman  signal  would  be  significantly 
weakened. 
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Table  1.  Band-gap  energies  in  units  of  eV  of  GaN,  AIN,  and  InN  [2], 


Structure 

GaN 

AIN 

InN 

Wurtzite 

3.39  (200  K) 

6.2  (300  K) 

1.89  (300  K) 

3.50  (1.6  K) 

6.28  (5  K) 

Zincblende 

3 .2-3.3  (300  K) 

5.11  (300  K) 

2.2  (300  K) 

3.  Raman  Scattering  of  GaN,  AIN,  and  InN  Films  and  Crystallites 

GaN,  AIN,  and  InN  based  materials  can  be  grown  in  the  wurtzite  as  well  as  the  zincblende 
structure,  depending  mainly  on  the  choice  of  substrates  and  growth  conditions  [10-12]. 
Unintentional  coexistence  of  both  structures  in  a  film  thus  may  occur.  This  section  focuses  on  the 
identification  of  the  structure  in  the  films  via  polarized  Raman  spectroscopy  and  analysis. 

3.1  Raman  Tensors  and  Structure  Identification  of  GaN,  AIN,  and  InN 

The  III-V  nitrides  are  highly  stable  in  the  hexagonal  WZ  structure  and  high  quality  material  in  that 
structure  has  been  achieved  [13].  However,  it  has  been  demonstrated  that  the  growth  of  films  in 
the  less  stable  cubic  ZB  structure  is  feasible  as  well  [14-17].  Due  to  the  polymorph  nature  of  the 
III-V  nitrides,  the  identification  of  material  structure  as  well  as  its  purity,  i.e.,  the  possibility  of 
polytype  material,  is  an  issue.  The  following  section  addresses  the  topic  of  Raman  spectroscopy 
for  structure  determination  of  GaN,  AIN,  and  InN  films  and  crystallites. 

The  WZ  crystal  structure  belongs  to  the  space  group  C64,. ,  and  the  group  theory  analysis 

predicts  the  zone-center  optical  modes  A\  +  2B\  +  E\  +  2Ei  [18].  The  A\,  E\,  and  the  two  Ej 
modes  are  Raman  active,  while  the  B\  modes  are  silent,  i.e.,  forbidden  in  Raman  scattering. 
Furthermore,  the  A\  and  E\  modes  are  polar:  their  vibrations  polarize  the  unit  cell  which  results  in 


5 


the  creation  of  a  long-range  electrostatic  field.  The  effect  of  this  field  manifests  itself  in  the 
splitting  of  the  A,,  E,  modes  into  longitudinal  optical  (LO)  and  transverse  optical  (TO) 
components,  thus  creating  the  A\(LO,TO)  and  E^LO.TO)  modes.  Figure  1  depicts  the  scheme  of 
the  vibrational  modes  in  the  wurtzite  structure  [19].  The  ZB  structure  belongs  to  the  space  group 
Tj  and  group  theory  predicts  one  Raman  active  mode  of  T2  representation:  it  is  a  polar  mode 

which  splits  into  the  TO  and  LO  components  [18,  20]. 

Tables  2  and  3  list  the  experimental  Raman  frequencies  that  have  been  reported  for  WZ-  and 
ZB-GaN,  respectively  [14,  16,  17,  21-28],  In  Table  4  the  Raman  frequencies  of  AIN  and  InN  are 
listed  as  well.  There  are  small  differences  in  the  Raman  frequencies  presented  in  Tables  2,  3  and  4, 
which  in  most  cases  are  due  to  differences  in  material  quality.  Moreover,  as  can  be  observed  in  the 
Tables,  the  Raman  frequency  of  the  TO  mode  of  the  ZB  lies  between  the  frequencies  of  the  E\ 
(TO)  and  the  A  i  (TO)  of  the  WZ;  similarly,  the  frequencies  of  the  LO ,  E\(LO),  and  the  A\(TO)  lie 
in  the  same  range.  The  close  proximity  of  the  ZB  and  the  WZ  phonons  has  been  attributed  to  the 
close  relation  between  these  two  structures  [29],  and  is  explained  as  follows.  The  WZ  structure, 
consisting  of  four  atoms  per  primitive  cell,  can  be  obtained  from  the  ZB  structure,  consisting  of 
two  atoms  per  cell,  by  a  rearrangement  of  the  atom-planes  that  are  perpendicular  to  the  (1 1 1)  axis 
[29].  Due  to  this  relation  the  WZ  phonon  dispersion  curve  along  [0001]  may  be  obtained  via 
folding  of  the  ZB  phonon  dispersion  along  the  [1 1 1]  direction,  thus  doubling  the  phonon  modes  in 
the  WZ  Brillouin  zone  as  required  by  the  increased  number  of  the  atoms  per  primitive  cell. 

Figure  2  presents  a  description  of  the  phonon  dispersion  of  the  two  structures.  In  the  figure, 
the  LO' s  and  the  TO' s  of  the  ZB  structure  are  plotted  along  the  T-L  [111]  direction  (thick  lines); 
the  folding  creates  the  two  E2  branches  and  the  two  B}  branches  (thin  lines)  along  the  T-A  [0001] 
direction  of  the  WZ  structure.  Moreover,  due  to  the  WZ  crystal  potential,  a  splitting  occurs  in 
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which  the  E\(TO)  and  A\(TO),  and  the  Ex(LO)  and  A\(LO)  branches  (dashed  lines)  are  created 
from  the  original  TO  and  LO  branches,  respectively  [29].  Thus,  as  can  be  seen  in  Fig.  2,  the 
folding  and  splitting  mechanisms  cause  the  Raman  frequencies  of  the  A](TO)  and  E\(TO)  modes 
of  the  WZ  structure  to  be  in  the  same  range  as  that  of  the  TO  mode  of  the  ZB  structure;  the  same 
holds  for  the  LO  modes. 


Table  2.  Raman  frequencies  (cm'1)  of  the  WZ-GaN  modes. 


E\ 

145 

144 

143 

144 

142 

E22 

568 

568 

569 

568 

570 

569 

567 

A\(TO) 

533 

533 

532 

533 

533 

530 

A  i  (LO) 

i 

738 

735 

E\(TO) 

559 

559 

560 

559 

561 

558 

E\(LO) 

726 

743 

Ref. 

[21] 

[22] 

[23] 

[25] 

[24] 

[26] 

[27] 

Table  3.  Raman  frequencies  (cm'1)  of  the  ZB-GaN  modes. 


TO 

554 

555 

555 

LO 

730 

740 

740 

737 

Ref. 

[16] 

[17] 

[28] 

[14] 
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Table  4.  The  Raman  frequencies  (cm'1)  of  WZ-A1N  and  WZ-  and  ZB-InN. 


WZ-A1N 

WZ-lnN 

ZB-A1N 

E\ 

241 

252 

246 

249 

655 (TO) 

E; 

660 

660 

655 

657 

495 

491 

Ai(TO) 

607 

614 

608 

610 

A\(LO) 

893 

890 

596 

590 

902  (LO) 

E\(TO) 

673 

668 

670 

E\(LO) 

924 

916 

913 

Ref. 

[30] 

[31] 

[27] 

[32] 

[33] 

[34] 

[15] 

The  identification  of  the  Raman  mode  may  be  achieved  using  the  Raman  tensors.  Equations 
(3)  and  (4)  present  the  Raman  tensors  for  the  ZB  and  WZ  structures  respectively  [18].  The  x, 
and  z  in  the  equations  represent  the  phonon-polarization  direction  for  each  of  the  polar  modes. 
Each  non-zero  entry  in  the  matrices  represents  a  coupling  of  the  polarization  of  the  incoming  and 
scattered  light  to  the  lattice  vibrations  that  results  in  an  allowed  Raman  spectral  line  via  the 
second  term  of  (1).  Throughout  this  review  Porto  notation  a(bc)d  will  be  used  to  describe  a 
Raman  geometry:  a  and  d  represent  the  wavevector  direction  of  the  incoming  and  scattered  light 
respectively,  while  b  and  c  represent  their  polarization  state  [35],  Additionally,  for  each  of  the 
polar  modes  the  direction  of  the  phonon  propagation  with  respect  to  its  polarization  direction 
determines  whether  the  observed  spectral  line  is  of  a  longitudinal  or  a  transverse  mode.  The 
phonon  propagation  direction  can  be  ascertained  via  the  wavevector  conservation  law:  k,.  =  ks  + 
q,  where  k,.  and  ks  are  the  wavevectors  of  the  incoming  and  the  scattered  light,  respectively,  and  q 
is  the  phonon  wavevector. 
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Thus,  for  example,  the  Raman  backscatteririg  configuration  z(xy)z  would  give  rise  to  the  LO 


mode  in  a  case  of  ZB-(OOl)  oriented  crystallographic  plane  or  to  the  two  E2  modes  in  a  case  of 


WZ-(OOOl)  plane.  The  other  modes  are  said  to  be  symmetry  forbidden  at  that  configuration  and 
wifi  not  appear  in  the  spectra  unless  a  mechanism  exists,  such  as  internal  reflection  due  to  defects 


and  interfaces,  that  may  cause  the  forbidden  scattering. 
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3.2  Wurtzite  and  Zincblende  Phases  of  GaN 

GaN  in  the  WZ  structure  has  been  grown  in  the  crystalline  form  of  needles  and  platelets  as  well  as 
heteroepitaxial  films  [1,  2,  5],  The  most  widely  used  substrates  are  sapphire(OOOl),  6H- 
SiC(OOOl),  and  ZnO(OOOl).  Their  thermal  expansion  coefficients  and  lattice  constants  are 
approximately  matched  to  GaN.  Among  the  first  to  identify  the  Raman  modes  in  WZ-GaN  were 
Manchon  et  al.  [21],  Lemos  et  al.  [22],  Bums  et  al.  [23]  and  Cingolani  et  al.  [25].  followed  by 
Murugkar  et  al.  [24]  and  Azuhata  et  al.  [26].  Table  5  lists  the  scattering  configurations  and  the 
observable  modes  for  the  WZ-GaN  film  studied  by  Azuhata  et  al.  [26].  The  2  pm  GaN  film  in  that 
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study  was  grown  on  3  sapphire(OOOl)  substrate  via  MOCVD  and  was  capped  with  a  thin  layer  of 
InGaN.  The  sample  quality  was  reported  to  be  comparable  to  that  of  a  double-hetero structure 
light-emitting  diode  [26].  Figures  3  and  4  show  the  Raman  spectra  of  that  film. 


Table  5.  Scattering  configurations  and 
observable  Raman  modes  of  WZ-GaN(OOOl)  [26]. 


Raman  geometry 

Raman  mode 

y(xr)y 

A]  (TO),  E\,  E; 

y{zz)y 

A\  (TO) 

y(zx)y 

E^c^o) 

z{xy)z 

El  El 

z(xx)5 

At  (LO),  El,  El 

The  values  of  the  LO-  and  the  719-phonon  may  be  used  to  determine  the  dielectric  constants 
of  a  material  via  the  Lyddane-Sachs-Teller  relation  [36] 


^Z  =  — 


(5) 


In  this  relation  So  is  the  static  dielectric  constant  of  the  crystal,  while  s.v  is  the  dielectric  constant 
at  optical  frequencies,  and  o\  and  ah  are  the  LO-  and  the  7<9-frequency,  respectively.  Since  GaN 


is  an  uniaxial  crystal,  there  are  two  values  for  each  of  the  dielectric  constants  defined  as  Sj  and  s„, 
respectively;  the  directions  given  are  relative  to  the  c-axis.  In  that  respect,  the  phonons  of  E\ 
symmetry  correspond  to  the  perpendicular  Lyddane-Sachs-Teller  relation  whereas  the  A  \  phonons 
to  the  parallel  relation.  Moreover,  infrared  reflectivity  measurements  have  indicated  that  for  GaN 
£j.cc  =  s„co,  i.e.,  in  the  optical  frequency  range  GaN  is  approximately  an  isotropic  material  [21]. 
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Using  the  Raman  frequencies  of  the  LO  and  TO,  Azuhata  et  al.  calculated  the  perpendicular  and 
the  parallel  static  dielectric  constants  of  the  GaN  film  to  be  £±o  -  9.28  and  e„o  =  10.1,  respectively 
(for  £_l~  =  £„:c  =  5.29)  [26].  These  values  are  in  accord  with  the  ones  obtained  via  infrared 
reflectivity  and  Kramers-Kronig  analyses  [37]. 

In  contrast  to  the  extensive  studies  that  have  been  conducted  on  WZ-GaN  material,  the  cubic 
phase  is  as  yet  a  new  and  largely  unexplored  area.  Raman  and  X-ray  scattering  studies  of  ZB-GaN 
films,  grown  on  GaAs(OOl)  substrates  by  metalorganic  vapor  phase  epitaxy  (MOVPE).  have  been 
reported  by  Miyoshi  et  al.  [16].  In  their  study,  the  structure  of  the  films  was  first  identified  via  X- 
ray  diffraction  and  then  confirmed  by  Raman  analysis.  The  X-ray  spectra  exhibit  the  (200) 
diffraction  peak  of  the  cubic  face  whereas  the  (0002)  diffraction  of  the  WZ  symmetry  is  absent.  In 
order  to  verify  the  X-ray  findings,  Raman  spectra  was  acquired  in  the  z(xy)z  configuration  for 
which  the  LO  of  the  cubic  symmetry  is  an  allowed  mode  from  the  (001)  face.  The  spectra  of  that 
study  consist  of  a  line  at  «  730  cm'1  which  has  been  attributed  by  Miyoshi  et  al.  to  the  LO  mode  of 
the  ZB-GaN  [16].  Moreover,  from  the  X-ray  data  the  authors  calculated  the  lattice  constant  of 
the  ZB-GaN  to  be  4.5  A,  which  agrees  well  with  their  theoretical  calculation.  In  light  of  that, 
Miyoshi  et  al.  concluded  that  despite  the  20%  lattice  mismatch  with  the  GaAs  substrate,  the  ZB- 
GaN  films  («  0.4  pm  thick)  are  almost  strain-free  due  to  a  relaxation  mechanism  occurring  at  the 
early  stages  of  the  growth  [16]. 

Phonon  characteristics  of  cubic,  hexagonal,  and  mixed  phase  GaN  films  grown  via  molecular 
beam  epitaxy  (MBE)  on  GaAs(001)  substrates  were  studied  by  Giehler  et  al.  [17].  Each  film  in 
this  study  was  grown  under  different  conditions  to  invoke  the  growth  of  the  different  structure 
types.  The  structural  properties  were  analyzed  by  high-energy  electron  diffraction  and  X-ray 
scattering,  and  the  phonon  characteristics  were  investigated  utilizing  infrared  transmission  and 
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Raman  spectroscopy.  The  structural  analysis  conducted  by  Giehler  et  al.  revealed  that  the  sample 
grown  under  Ga-rich  conditions  (sample  A)  consists  of  nonepitaxial  hexagonal  columns  whose  c- 
axes  are  aligned  parallel  to  the  [001]  direction  of  the  substrate.  The  film  grown  under  N-rich 
conditions  (sample  B)  exhibits  a  phase  mixture,  while  the  film  grown  under  near-stoichiometric 
conditions  (sample  C)  was  mainly  cubic. 

The  infrared  transmission  spectra  of  sample  A  consist  of  signals  at  558  cm 1  and  735  cm  , 
which  have  been  attributed  by  Giehler  et  al.  to  the  Ex(TO)  and  Ax(LO)  phonons,  respectively,  of 
the  WZ-GaN.  The  spectra  of  sample  C  reveals  transmittance  signals  at  552  cm 1  and  739  cm  that 
were  assigned  to  the  TO  and  LO  phonons,  respectively,  of  the  ZB-GaN  [17].  The  spectra  in  the 
TO  frequency  range  of  the  mixed  sample  exhibit  a  broad  signal  centered  approximately  between 
the  E](TOj  and  the  TO  modes;  similarly,  the  signal  in  the  LO  range  is  broad  and  centered  between 
the  A\(LO)  and  the  LO  modes.  The  Raman  spectra  of  the  three  samples  are  presented  in  Fig.  5; 
the  spectra  were  acquired  in  the  z{xy)z  configuration.  The  spectra  of  sample  A  exhibit  the  E2 
mode  at  569  cm'1  and  a  broad  band  ~  680  cm'1  which  has  been  assigned  by  the  authors  to  a 
disordered  state.  The  spectra  of  sample  C  consist  of  the  allowed  LO  mode  at  740  cm'1  and  the 
symmetry  forbidden  TO  mode  at  554  cm'1.  The  forbidden  scattering  has  been  explained  by  Giehler 
et  al.  to  occur  via  multiple  reflections  inside  the  film.  The  spectra  of  the  mixed  phase  film  exhibit 
the  E]  of  the  WZ  structure,  the  LO  and  the  TO  of  the  ZB,  and  a  broad  band  due  to  a  disorder- 

activated  scattering  [17]. 

Comparative  Raman  studies  of  cubic  and  hexagonal  GaN  have  been  reported  also  by  Tabata 
et  al  [14].  In  that  study,  both  WZ  and  ZB  films  were  grown  on  GaAs(OOl)  substrates  via  MBE. 
The  nucleation  of  the  two  structures  was  achieved  by  varying  the  Ga-to-N  flux  ratio  in  the  initial 
stage  of  the  growth:  the  growth  of  the  ZB  film  was  initiated  using  N-rich  conditions,  and  that  of 
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the  WZ  phase  was  obtained  under  a  slight  excess  of  Ga  [14].  The  X-ray  spectra  of  the  sample 
grown  under  Ga-rich  condition  (sample  A)  exhibit  the  (0002)  and  the  (0004)  diffraction  peaks 
associated  with  the  WZ-GaN,  whereas  the  spectra  of  the  sample  grown  under  N-rich  condition 
(sample  B)  exhibit  the  (002)  peak  of  the  ZB-GaN.  From  these  findings,  the  authors  concluded 
that  sample  A  has  a  WZ  structure  with  the  surface-normal  parallel  to  the  c-axis,  and  sample  B 
consists  of  ZB  structure  with  the  surface-normal  parallel  to  the  cubic  axis  [14]. 

The  Raman  spectra  of  samples  A  and  B  are  presented  in  Figs.  6  and  7,  respectively.  The 
spectra  were  acquired  in  nearly  backscattering  geometry  with  polarized  incident  light  and 
unanalyzed  scattered  light.  Moreover,  the  spectra  in  Tabata's  study  were  acquired  at  three 
different  incident  wavelengths  in  order  to  enhance  the  Raman  signal  [14].  Figure  6  exhibits  the 
allowed  E\  mode  at  571  cm’1  as  well  as  the  Ax(LO)  mode  at  737  cm’1,  which  is  also  allowed  in 
that  configuration.  The  Raman  spectra  of  the  ZB-GaN  in  Fig.  7  exhibit  the  allowed  LO  mode  at 
741  cm'1,  and  the  forbidden  TO  mode  appears  as  well  at  555  cm'1.  The  appearance  of  the  TO 
mode  was  attributed  by  Tabata  et  al.  to  the  presence  of  a  short-range  perturbation  in  the  film, 
which  relaxes  the  Raman  selection  rules  [14], 

3.3  Wurtzite  and  Zincblende  Structure  of  AIN  and  InN 

In  the  following,  a  review  of  Raman  scattering  studies  in  AIN  and  InN  is  presented.  AIN  and  InN, 
like  GaN,  in  the  stable  state  have  a  wurtzite  crystal  structure  and  belong  to  the  space  group  C64v . 
The  cubic  zincblende  phase  of  AIN  and  InN  can  be  grown  as  well;  however,  these  materials  have 
been  studied  much  less  than  the  wurtzite  polytype.  The  Raman  selection  rules  leading  to  mode 
assignments  and  polytype  identification  are  the  same  as  those  for  GaN,  as  described  in  Section 

3.1. 
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Figure  8  shows  typical  Raman  spectra  of  'a  WZ-A1N  crystallite  grown  via  the  sublimation 
method  [27].  The  spectra  were  acquired  in  a  backscattering  geometry  from  the  a- face  of  the 
crystallite,  and  the  Raman  frequencies  are  listed  in  Table  4.  Figures  9  and  10  present  the  Raman 
spectra  of  a  ZB-AJN  film  along  with  that  of  ZB-AlGaN  films;  the  symmetry-allowed  LO  and  the 
forbidden  TO  modes  of  the  cubic  ZB  are  both  present  in  the  spectra  [15].  These  cubic  films  were 
grown  by  Harima  et  al.  via  the  MBE  method  on  cubic  SiC/Si  substrates,  and  the  cubic  phase  was 
confirmed  by  X-ray  diffraction  as  well  to  be  the  major  structure  [15]. 

Due  to  the  difficulties  in  achieving  reasonably  good  quality  InN  films,  information  on  its 
material  properties  is  as  yet  scarce.  Raman  scattering  studies  of  InN  films  grown  on  (0001) 
sapphire  substrates  via  metalorganic  vapor  phase  epitaxy  were  first  reported  by  Kwon  et  al.  [33]. 
X-ray  diffraction  indicated  that  the  structure  of  the  InN  films  was  wurtzite.  Identification  of  the 
Raman  modes  was  carried  out  via  the  Raman  selection  rules  in  a  backscattering  geometry;  the 
modes  are  listed  in  Table  4.  The  authors  reported  a  linewidth  of  20  cm'1  for  the  E]  mode,  a  value 
that  is  larger  compared  to  the  E]  linewidths  of  WZ-GaN  and  AIN,  which  usually  range  from  3  to 
8  cm'1.  Lee  et  al.  investigated  the  unusual  line  broadening  of  the  InN  Raman  modes  as  a  function 
of  growth  temperature  [34].  These  films  were  grown  at  temperature  ranges  of  325-600°C  via  the 
same  method  used  by  Kwon  et  al.  It  was  found  that  in  the  elevated  temperature  range  (450°C  and 
above),  the  E]  linewidth  is  *  6  cm'1  while  in  the  lower  range  the  linewidth  is  ~  70  cm'1;  the 
broadening  of  the  spectral  lines  was  accompanied  by  a  frequency  shift  as  well.  Lee  et  al.  suggested 
that  this  significant  line  broadening  reflects  the  coexistent  of  cubic  and  amorphous  phases  in  the 
WZ  structure  of  the  InN  films  [34], 
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4.  Stress  Analysis  and  Substrate  Issues  for  Epitaxial  Growth 

Among  the  important  factors  leading  to  the  synthesis  of  high  quality  films  is  the  availability  of 
matching  substrates.  In  order  to  minimize  the  internal  stress,  the  substrate  under  consideration 
should  have  lattice  parameters  and  thermal  expansion  coefficients  similar  to  these  of  the  film. 
Table  6  lists  some  of  the  substrates  most  commonly  used  in  the  growth  of  III-V  nitride  films,  as 
well  as  their  lattice  parameters  and  thermal  expansion  coefficients  [2]. 

One  of  the  most  informative  methods  of  measuring  stress  in  materials  is  Raman  spectroscopy 
[19,  38,  39];  however,  it  requires  an  a-priori  knowledge  of  the  Raman  pressure  coefficients  that 
relate  the  peak  position  to  the  stress.  The  following  reviews  some  of  the  Raman  studies 
concerning  the  determination  of  the  pressure  coefficients  and  the  stress  state  in  GaN  and  AIN. 


Table  6.  Material  properties  for  the  GaN-and  AlN-substrate  system  [2]. 


Material 

Lattice  parameters  (A) 

Coefficients  of  thermal 
expansion  ( 1  O^/K) 

6H-SiC 

o  =  3.080 

4.2 

c  =  15.12 

4.68 

Sapphire 

a  =  4.758 

7.5 

c  =  12.99 

8.5 

ZnO 

0  =  3.252 

2.9 

c  =  5.213 

4.75 

WZ-AIN 

0  =  3.112 

4.2 

c  =  4.982 

5.3 

WZ-GaN 

o  =  3.189 

5.59 

c  =  5.185 

3.17 

4.1  Stress  Analysis  of  GaN  Films 

A  detailed  study  of  the  effect  of  pressure  on  the  crystalline  WZ-GaN  Raman  spectral  lines  was 
reported  by  Perlin  et  al.  [38].  The  applied  pressure  was  hydrostatic  in  the  range  0  to  50  GPa.  The 
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pressure  dependence  of  the  Raman  frequencies  of  the  A\(TO),  E\(TO),  E{ ,  and  E\  modes  was 
shown  to  follow  the  quadratic  relation  [38]: 

n^cm'1)  =  ftb  +  G\P  +  cf>P~  (^) 

where  P  is  the  applied  pressure  in  units  of  GPa,  ccc,  is  the  Raman  frequency  at  zero  applied 
pressure,  and  cr,  and  05  are  the  first  and  second  order  pressure  coefficients,  respectively.  The 
hydrostatic  pressure  coefficients  derived  from  first  principal  density-functional  calculations  by 
Gorczyca  et  al.  [19]  were  found  to  agree  well  with  the  experimental  results  [38],  Table  7 
summarizes  the  set  of  parameters  for  the  various  Raman  modes. 


Table  7.  The  hydrostatic  pressure  coefficients  of  WZ-GaN,  o 1  and  02,  in  units  of  cm  '/GPa  and 
cm'V(GPa)2,  respectively,  as  obtained  from  theory  [19]  and  experiment  [38]. 


Raman  mode 

Theory 

Experiment 

Cl 

02 

Cl 

Oi 

E\(TO) 

4.10 

0.013 

3.68 

-0.0078 

E\ 

-0.15 

0.006 

-0.25 

-0.0017 

E; 

4.46 

0.018 

4.17 

-0.0136 

A, (TO) 

4.08 

0.024 

4.06 

-0.0127 

The  stress-strain  relation  for  an  isotropic  material  under  hydrostatic  pressure  is  [40,  41]: 

s  =  F]cr(\-2v),  C7-1) 

where  sis  the  strain,  E  the  Young’s  modulus,  crthe  stress,  and  v  the  Poisson  ratio.  However,  for 
thin  epitaxial  films  the  stress  is  more  likely  to  be  two  dimensional  due  to  film-substrate  lattice  and 
thermal  mismatches  (see  Table  6).  For  the  case  of  biaxial  stress,  <ja,  a  strain  is  induced  in  the  basal 
plane,  sa,  and  along  the  c-direction,  sc;  the  isotropic  constituent  equations  are  [41]: 
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(7.2a) 


ft=  E]aa(\-2v) 

sc  =  -K]aa2v.  (7-2b) 

For  thin  epitaxial  films,  the  biaxial  stress  approximations  in  (7.2  a,  b)  may  be  used  to  describe  the 
stress  state  of  a  thin  film  if  the  stress  contribution  from  point  defects  is  negligible  [40]. 

Kisielowski  et  al.  have  analyzed  the  stress  state  of  GaN  films  grown  on  the  c-face  of  SiC  and 
sapphire  [40].  The  authors  analyzed  the  stress  in  the  films  in  terms  of  the  biaxial  isotropic 
approach  described  in  (7.2).  It  was  shown  that,  due  to  point  defects,  an  additional  hydrostatic 
component  has  to  be  superimposed  to  the  biaxial  component  in  order  to  fully  describe  the  stress 
state  of  the  films.  Moreover,  Kisielowski  et  al.  determined  that  a  biaxial  stress  of  1  GPa  would 
shift  the  Raman  peak  position  of  the  E \  mode  by  «  4.2  cm',  a  result  consistent  with  the  values 

listed  in  Table  7. 

Strain  effect  studies  in  epitaxial  GaN  grown  on  AlN-buffered  Si(lll)  have  been  reported  by 
Meng  et  al.  [42].  One  of  the  objectives  of  the  study  was  the  determination  of  the  Raman  stress 
coefficient  for  the  biaxial  stress  state,  taking  into  account  the  hexagonal  anisotropy.  An  outline  of 
their  method  is  presented  here.  The  films  were  grown  via  RF-glow  discharge  reactive  magnetron 
sputtering  at  various  RF  input  powers.  Raman  and  X-ray  spectroscopy  established  that  the  film 
structure  is  WZ  with  an  (0001)  orientation.  In  their  analysis,  Meng  et  al.  assumed  a  biaxial  stress 
state  and  took  into  consideration  the  hexagonal  anisotropy;  the  constituent  relation  in  that  case  is 


Using  the  elastic  constants  Q  for  WZ-GaN  in  (8)  the  authors  obtained  aa  =  -202 sc.  The  values  of 
the  vertical  strain  sc  determined  via  X-ray  measurements  were  correlated  to  the  RF  input  power 


17 


and  to  the  E]  Raman  peak  position;  in  both  cases  a  direct  correlation  was  found.  The  vertical 
strain-Raman  frequency  relation  found  for  the  E2  mode  was  o/sc)  —  561+701^,  which  in 
conjunction  with  (8)  yields  the  relation  between  the  biaxial  stress  and  the  Raman  frequency  to  be 
7a)  =  561+4. 4cra.  Thus  the  biaxial  pressure  coefficient  found  in  the  study  of  Meng  et  al.  (4.4 
cnf'/GPa)  [42]  is  comparable  to  the  hydrostatic  pressure  coefficient  [38]. 

Other  studies  of  stress-related  phenomena  have  been  reported  [43-45];  a  brief  summary  is 
given  below.  Thermal  stress  effects  on  GaN  films  of  differing  thickness  grown  on  the  c-plane  of 
sapphire  substrates  were  investigated  by  Kozawa  et  al  [45].  As  can  be  seen  in  Table  6,  a 
significant  difference  exists  between  the  thermal  expansion  coefficients  of  GaN  and  those  of 
sapphire,  which  produces  a  compressive  stress  in  the  film.  In  the  studies  of  Kozawa  et  al.,  the 
magnitude  of  the  biaxial  compressive  stress  in  each  of  the  GaN  films  was  obtained  via  the 
curvature  wafer  bending  method;  the  obtained  values  were  correlated  to  the  peak  positions  of  the 
E l  Raman  mode.  It  was  found  that  the  Raman  pressure  coefficient  is  6.2  cm'VGPa,  a  value 
somewhat  greater  than  the  one  listed  in  Table  7.  Additionally,  the  authors  demonstrated  that  the 
compressive  stress  decreases  as  the  thickness  of  the  film  increases:  for  ~  5  pm  and  50  pm  films 
the  Raman  peak  position  was  found  to  be  at  569  and  567.5  cm’1,  respectively  [45]. 

Rieger  et  al.  investigated  the  influence  of  the  AIN  buffer  layer,  when  deposited  on  the  c-plane 
of  sapphire  substrates,  on  the  stress  state  of  the  GaN  films  [43].  In  order  to  estimate  the  stress 
magnitude,  they  analyzed  the  Raman  peak  position  of  the  E2  mode  as  a  function  of  the  AIN  layer 
thickness.  Rieger  et  al.  observed  a  pronounced  reduction  of  the  compressive  stress  in  the  GaN 
films  with  increasing  ALN  buffer  layer  thickness.  According  to  the  authors,  an  AIN  layer  thicker 
than  200  nm  eliminates  the  compressive  stress  completely  [43],  Depth  profiles  of  the  strain  in  a 
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220-pm  thick  GaN  film  grown  on  sapphire  were  investigated  by  Siegle  et  al.  utilizing  Raman 
spectroscopy  PL  and  CL  imaging  [44],  The  depth  profile  of  the  stress  was  determined  from  the 
peak  position  of  the  E \  mode.  The  authors  found  that  the  stress  decreases  nearly  exponentially 
with  increasing  distance  from  the  substrate;  the  film  was  found  to  be  relaxed  at  a  distance  ~  35 
pm  away  from  the  substrate. 

4.2  Stress  Analysis  in  WZ-A1N 

The  influence  of  hydrostatic  pressure  on  the  Raman  frequencies  of  WZ-A1N  has  been  investigated 
by  various  groups  [19,  30,  46].  The  Raman  pressure  coefficients  for  the  AIN  modes,  defined  in 
(7.1),  are  listed  in  Table  8. 

Table  8.  The  hydrostatic  pressure  coefficients  of  WZ-A1N,  cr,  and  ch.  in  units  of  cm'VGPa  and 
cm'V(GPa)2  respectively  as  obtained  from  theory  [19]  and  experiment  [30]. 


Theory 

Experiment 

Raman  mode 

V\ 

02 

cr, 

Ch 

Ei(TO) 

4.36 

0.059 

E\ 

-0.29 

0.022 

El 

4.79 

0.063 

3.99 

0.035 

A\(TO) 

4.29 

0.019 

4.63 

-0.01 

E\(LO) 

1.67 

0.27 

Figures  1 1  and  12  [19]  show  the  functional  behavior  of  the  Raman  modes  under  pressure;  the 
dots  represent  the  experimental  data  of  [30]  while  the  lines  represent  the  model  calculation  [19]. 
One  of  the  most  intriguing  aspects  of  these  results  is  the  anomalous  behavior  of  the  E\  mode  of 
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AIN  as  well  as  of  GaN;  unlike  the  other  modes  it  decreases  to  lower  frequencies  with  increasing 
pressure.  This  tendency  is  better  expressed  in  terms  of  the  Gruneisen  parameter  y,  [47]: 


r,= 


B0  (  d<y, 


co, 


V  d  P 


(9) 


where  the  s  are  the  given  phonon  modes,  B0  is  the  bulk  modulus,  and  P  is  the  pressure.  Thus  all 
the  AIN  and  GaN  modes  have  positive  y  values,  excepting  the  E\  mode,  which  has  a  negative  y 
value.  Similar  behavior  has  been  found  in  CdS  and  ZnO,  both  of  which  have  a  WZ  structure  [47]. 
It  has  been  suggested  that  the  negative  y  value  may  indicate  a  softening  of  the  lattice  to  a 
particular  mode  under  compression,  preceding  a  first-order  transition  [19,  47].  Moreover,  it  has 
been  established  that  WZ-ZnO  and  -CdS  undergo  a  phase  transition  to  the  NaCl  structure  [47]; 
the  same  phase  transition  has  been  suggested  to  occur  in  GaN  crystallites  at  *  45  GPa  [38]. 
However,  the  reason  for  the  softening  of  the  E\  mode  in  the  wurtzite  structure  materials  needs 
further  investigation. 


5.  Raman  Analysis  of  the  Quasi-Modes  in  AIN 

The  theory  developed  by  Loudon  formulates  that  in  uniaxial  materials  the  polar  phonon 
characteristics  may  be  affected  via  two  interaction  mechanisms:  one  due  to  the  long  range 
electrostatic  field,  and  the  other  due  to  the  short  range  field  which  exhibits  the  anisotropy  of  the 
vibrational  force  constants  [18,  48].  The  phonon  dynamics  and  thus  the  Raman  spectra  depend  on 
which  of  the  two  mechanisms  is  the  dominant  interaction.  Figure  13  depicts  the  Raman  frequency 
scheme  and  the  direction  of  phonon  vibrations  for  both  interactions  [49]. 

For  the  case  where  the  long  range  electrostatic  field  is  the  dominant  mechanism,  the 
interaction  of  the  polar  phonons  with  the  long  range  electrostatic  field  may  result  in  a  significant 
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frequency  separation  between  the  group  of  the  TO  phonons  relative  to  that  of  the  LO  phonons. 
Moreover,  the  TO  phonons  belonging  to  different  symmetry  are  grouped  together  in  a  relatively 
narrow  frequency  range;  the  same  holds  for  the  LO  phonons.  One  consequence  of  the  dominant 
electrostatic  field  interaction  is  that  under  certain  propagation  and  polarization  conditions, 
phonons  of  mixed  A\  and  E\  symmetry  character  exist  and  can  be  observed  in  the  Raman  spectra 
[48].  These  mixed  symmetry  modes  are  termed  quasi-TO  and  -TO  modes.  The  frequencies  of 
such  quasi  modes  are  predicted  by  the  theory  to  be  between  the  values  of  the  pure  A,  and  E\ 
modes  for  each  of  the  LO  and  the  TO  bands.  Alternatively,  for  the  case  where  the  short  range 
interatomic  forces  are  dominant,  the  LO-TO  splitting  will  be  small  and  in  this  case  the  TO  and  LO 
mode  of  each  symmetry  group  will  occur  in  a  relatively  narrow  band. 

In  AIN  the  E{(TO)  and  A\(TO)  Raman  frequencies  are  grouped  together  within  about  60  cm'1, 
the  E]  (LO)  and  A ,  (LO)  within  20  cm'1,  and  the  LO-TO  group-splitting  is  «  220  cm'1.  The 
frequency  scheme  of  AIN  thus  implies  the  dominance  of  the  long-range  electrostatic  force 
interaction.  The  mode-mixing  in  AIN  may  occur  if  the  propagation  direction  (^-vector)  of  the 
quasi-polar  phonons  lies  in  the  plane  which  spans  between  the  c  and  the  a\  (or  a2)  crystallographic 
axes  [50].  This  plane  is  referred  as  to  the  mixing  plane.  For  example,  the  pure  A\  phonon  has  a  c- 
direction  polarization  while  the  pure  E\  phonon  is  polarized  in  the  basal  plane;  thus  a  quasi¬ 
phonon  with  the  ^-vector  between  the  c  and  the  ax  axes  would  exhibit  a  mixed  polarization  of  Ar 
E\  symmetry.  When  the  q-ve ctor  lies  along  the  crystallographic  axes  or  in  the  a\-aj  plane,  only 
pure  phonons  are  observed  in  the  spectra.  A  detailed  analysis  of  the  Raman  configurations  that 
enable  the  observation  of  the  quasi  modes  in  WZ  materials  is  given  in  [50]  and  [49]. 
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Raman  studies  of  quasi  modes  in  AIN  have  been  reported  by  Filippidis  et  al.  [32]  and  Bergman 
et  al.  [49].  In  these  studies  the  ^-vector  direction  was  achieved  by  rotating  the  crystallite  by  an 
angle  <9  relative  to  a  fixed  laboratory  coordinate  system,  as  depicted  in  Fig.  14.  Figures  15  and  16 
present  Raman  spectra  for  various  rotations  in  a  scattering  configuration  that  enables  the  quasi- 
TO  and  -LO  to  be  observable  [49].  Figure  17  shows  the  behavior  of  the  quasi-modes  as  a 
function  of  the  phonon  propagation  angle  J3,  where  /?  is  related  to  6  via  Snell’s  law  and  a  small 
anisotropy  is  assumed.  In  Fig.  17  the  points  represent  the  experimental  results  while  the  lines 
represent  Loudon’s  model  [48,  50]: 

CDqvo,  =  Co\ia0)  cos 2  (90  -  p)+co;4im  sin- (90  -  p)  (10a) 

co2Q(m  =co2A,(LO)cos2(90- p)+co2El(U)}sin2(90- p) .  (10b) . 

As  was  discussed  in  [49].  one  important  implication  of  these  results  involves  the  accurate 
determination  of  the  Raman  frequency  acquired  from  AJN  thin  films.  In  such  films,  due  to  their 
pm-size  dimension  along  one  of  the  crystallographic  axes,  the  determination  of  the  various  Raman 
frequencies  requires  the  data  to  be  acquired  at  grazing  angle  geometry.  Such  experimental  setup, 
which  deviates  by  an  angle  from  the  crystallographic  coordinate  system,  may  result  in  a  mode 
mixing  of  a  characteristic  Raman  frequency.  Due  to  the  strong  angular  dispersion  of  the  quasi -TO 
mode  frequency  (see  Fig.  17),  a  disparity  in  the  TO  Raman  frequency  values  of  thin  films  is 
expected  to  be  observed;  less  disparity  is  expected  for  the  LO  mode.  A  similar  effect  in  which  the 
frequency  is  shifted  from  its  pure  value  may  also  occur  in  films  grown  at  off-crystallographic 
directions.  In  measuring  the  stress  state  of  the  films,  the  possible  additional  shift  due  to  the 
electrostatic  coupling  needs  to  be  considered. 
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6.  Phonon-Plasmon  Interaction  in  GaN  Films  and  Crystallites 

Raman  spectroscopy  has  proven  to  be  a  useful  tool  in  analyzing  the  effect  of  free  electrons  on  the 
lattice  dynamics  of  GaN  which,  as  grown,  tends  to  be  an  w-type  material  [4],  A  free  carrier 
concentration  in  GaN  may  be  achieved  via  intentional  n-type  doping  or  via  the  growth  method  [5, 
51].  In  general,  when  an  appreciable  free  carrier  concentration  is  present  in  a  polar- 
semiconductor,  a  coupling  of  the  IO-phonons  to  the  plasma  oscillations  of  the  free  carriers 
(plasmons)  may  occur  [52],  The  phonon-plasmon  interaction  results  in  a  characteristic  Raman 
scattering  that  may  yield  information  on  the  free  carrier  density  of  a  given  sample  [51,  53-58]. 
This  section  examines  the  topic  of  Raman  phonon-plasmon  coupled  modes  in  GaN. 

The  determination  of  the  free  carrier  concentration  via  Raman  spectroscopy  is  often  more 
advantageous  than  that  via  electrical  measurements  since  no  contacts  are  required.  The  plasma 
frequency  CQp  at  wavevector  a  =  0  is  related  to  the  free  carrier  concentration  n  through  the  relation 
[36]: 

a2  Anne2/  (11) 

r  /™ 

where  £<*  is  the  high  frequency  dielectric  constant,  e  is  the  charge,  and  m*  is  the  effective  mass  of 
the  free  carriers.  Moreover,  in  the  case  when  the  plasma  oscillations  are  not  overdamped  the 
frequencies  of  the  coupled  plasmon-IO-phonon  may  be  approximated  by  [59]: 

°>1  -4 <°WnY}  ■  <12) 

Therefore  the  Raman  spectra  in  principle  should  exhibit  two  spectral  lines  corresponding  to  an 
upper  branch  coupling  co-  and  a  lower  branch  co.  (also  known  as  L t  and  I  )  whose  frequencies 
depend  on  the  free  carrier  concentration.  No  10-mode  is  expected  in  the  spectra.  The  two 
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coupled  frequencies  have  been  observed  in  the  Raman  spectra  of  n-type  GaAs  as  well  as  other  III- 
V  semiconductors  of  the  ZB  structure  [20]. 

In  the  case  of  overdamped  plasmons  (12)  is  not  applicable;  the  spectra  cannot  be  analyzed  in 
terms  of  the  co,  and  ox  coupled  modes.  However,  the  case  of  large  damping  of  the  plasma 
oscillations  in  6H-SiC  (doped  with  nitrogen  ~  1019  cm'3)  has  been  analyzed  and  modeled  by  Klein 
et  al.  [60].  In  their  study  the  Raman  spectra  of  the  SiC  do  not  exhibit  the  ax,  and  ax  branches. 
Instead,  the  influence  of  the  coupling  manifests  itself  in  the  lineshape  of  the  ZO-mode:  it  is 
asymmetrically  broadened  and  slightly  shifted  toward  the  frequency  range  where  the  axr  was 
supposed  to  be  (in  the  case  of  small  damping).  Moreover,  Klein  et  al.  demonstrated  that  the  line 
broadening  as  well  as  intensity  reduction  is  directly  correlated  to  the  increase  of  carrier 
concentrations.  The  influence  of  damping  effect  on  the  Raman  scattering  efficiency  of  the  coupled 
modes  has  also  been  investigated  by  Hon  et  al.  [61]  and  by  Irmer  et  al.  [62]  in  the  case  of  n-type 
GaP.  The  Raman  spectra  were  found  to  exhibit  characteristics  similar  to  those  of  n-type  SiC.  The 
model  developed  to  explain  the  effect  of  the  large  damped  plasmons  on  the  coupled  modes  can  be 
expressed  as  [60-62] : 


s(ox)  =  e(<x>) 


1  + 


I(o))  *  [\,C,C2]Im(-\/e) 

col-co]  <»\ 

ox2  -  ox2  -  iyox  ox2  +  iypox 


(13a) 

(13b) 


where  I(co)  is  the  Raman  cross  section,  s(ox)  is  the  dielectric  function,  and  [l.C.C2]  is  the 
interference  factor  expressed  in  term  of  the  Faust-Henry  coefficient  C  [61].  In  (13b)  oi,  ox,  and  C0p 
are  the  frequencies  of  the  longitudinal  phonon,  the  transverse  phonon,  and  the  plasmons, 
respectively.  Additionally,  y  and  yp  are  the  damping  constants  of  the  phonons  and  the  plasmons, 
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respectively.  Equation  (13a),  when  fully  expanded,  predicts  the  Raman  lineshape  of  the  phonon- 
plasmon  coupled  modes. 

In  a  manner  similar  to  SiC  the  plasma  oscillations  in  GaN  are  considered  to  be  overdamped. 
Kozawa  et  al.  investigated  the  A\(LO)  phono n-plasmon  coupled  mode  in  WZ-GaN  film  as  a 
function  of  a  relatively  low  level  of  Si-dopant  concentration  [51].  The  dopant  concentration  in 
these  samples  ranged  from  about  lxlO16  to  2xl018  cm'3.  Figure  18  depicts  the  Raman  spectra  of 
the  A](LO)  mode  at  several  carrier  densities;  here  the  Raman  band  shifts  towards  the  high 
frequency  side  as  well  as  broadens  and  weakens  as  the  carrier  concentration  increases.  The 
spectra  of  the  A\(LO)  coupled  mode  as  well  as  the  model  predicted  by  the  lineshape  (13)  are 
presented  in  Fig.  19.  The  inset  to  the  figure  lists  the  values  of  the  fitting  parameters,  among  which 
are  the  plasma  damping  constant  y=  400  cm'1  and  the  plasma  frequency  Op  =1 19  cm'1,  the  latter 
from  which  the  carrier  concentration  may  be  evaluated  via  (11).  Figure  20  shows  the  carrier 
concentration  n  calculated  from  the  Raman  data  versus  the  one  obtained  from  the  Hall 
measurements.  As  shown  in  the  figure,  the  values  of  the  concentrations  obtained  in  both 
techniques  agree  [51]. 

Kirillov  et  al.  studied  the  effect  of  free  carriers  on  phonon-plasmon  interactions  in  WZ-GaN 
films  in  the  high  concentration  regime  [54],  The  concentrations  in  their  study  were  determined  via 
Hall  measurements  to  be  4x1 019  cm'3  for  sample  A  and  8xl019  cm'3  for  sample  B.  From  those 
measurements  the  plasma  damping  constant  for  sample  A  was  been  determined  to  be  1116  cm'1. 
Figure  21  presents  the  Raman  spectra  acquired  in  a  backscattering  geometry  from  the  c- face  of 
the  two  films:  neither  spectra  exhibits  the  allowed  A\(LO)  mode  at  «  735  cm'1.  In  order  to  gain 
further  insight  into  the  dynamics  of  the  phonon-free  carrier  interactions,  Kirillov  et  al.  calculated 
the  expected  spectra  using  Klein’s  model  expressed  in  (13a,b).  Figure  22  shows  the  calculated 
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spectra  for  a  free  carrier  concentration  of  4xl019  cm'3  for  two  values  of  the  plasma  damping 
constant.  The  dotted  line  in  Fig.  22  is  a  calculated  spectrum  for  a  hypothetical  sample  with  a  low 
damping  constant  of  100  cm'1.  This  value  is  typical  for  plasmons  in  GaAs  films;  in  that  material 
the  upper  branch  coupling  co-  and  a  lower  branch  ox  have  been  clearly  identified.  The  solid  line  in 
Fig.  22  is  the  calculated  spectrum  where  the  damping  constant  is  taken  to  be  that  of  sample  A 
(1116  cm'1).  As  shown  in  the  figure  the  lines  become  much  broader  as  the  plasma  damping 
constant  increases,  with  the  upper  branch  being  the  most  affected. 

In  contrast  to  the  study  of  Kozawa  et  al.  in  the  low  dopant  regime,  the  study  of  Kirillov  et  al. 
indicates  that  for  high  free  carrier  concentrations  Raman  spectroscopy  may  not  be  a  useful  tool  to 
extract  qualitative  information  on  the  concentration  [54],  This  is  because  the  line  of  the  upper 
branch  of  the  phonon-plasmon  coupling  is  too  broad  to  allow  a  meaningful  analysis.  However,  the 
model  calculations  presented  in  Fig.  22  reveal  that  the  line  of  the  lower  branch  at  *  530  cm'1  is  not 
too  broad  and  thus  potentially  may  be  used  in  the  concentration  analysis.  This  line  corresponds  to 
the  523  and  530  cm'1  spectral  lines  of  samples  A  and  B,  respectively  (see  Fig.  21).  However,  due 
to  defects  and  inhomogeneities  in  the  film,  which  relax  the  Raman  selection  rules,  the  A\(TO) 
mode  near  530  cm'1  may  be  present  in  the  spectra  as  well.  Thus  the  identification  of  the  lower 
phonon-plasmon  coupling  mode  has  to  be  dealt  with  cautiously. 

A  detailed  investigation  into  the  issue  of  the  lower  branch  of  the  phonon-plasmon  coupled 
modes  in  the  case  of  high  carrier  concentration  was  carried  out  by  Demangeot  et  al.  [55].  The 
GaN  sample  was  found  via  infrared  reflectivity  measurements  to  contain  a  high  level  of  free 
carrier  «  102°  cm'3,  and  the  plasma  damping  constant  was  determined  to  be  1480  cm'1  [55],  Figure 
23  shows  Raman  spectra  of  the  GaN  film  that  were  acquired  at  two  locations  on  the  sample.  As 
can  be  seen  in  the  figure,  a  relatively  broad  Raman  peak  at  525  cm 1  is  present  in  the  spectra. 
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Demangeot  et  al.  used  the  model  of  (13a,  b)  and  found  a  good  fit  (dotted  line  in  the  figure)  to  the 
Raman  data  of  the  low  branch  coupled  mode.  As  can  be  observed  in  Fig.  23,  the  Raman  spectra  in 
the  high  frequency  range  do  not  exhibit  the  A{(LO)  mode  nor  a  pronounced  high-branch  coupled 
mode,  results  which  are  consistent  with  the  model  calculation.  Additionally,  the  results  of 
Demangeot  et  al.  indicate  that  the  distribution  of  the  free  carriers  is  not  uniform  across  the  GaN 
sample. 

The  fact  that  in  the  limit  of  low  free  carrier  concentration  the  modified  A\(LO)  is  present  in  the 
Raman  spectra  led  Wetzel  et  al.  to  quantitatively  determine  the  concentrations  [57].  In  this  study 
they  investigated  bulk  GaN  crystallites  and  correlated  the  free  carrier  concentration  to  the  Raman 
peak  position  of  the  coupled  T,(Z,0)-plasmon  mode.  The  authors  suggested  that  the  following 
correlation  holds  in  the  concentration  range  of  1017  cm J  <n  <  10  9  cm  [57]: 

n(cd)  =  l.lxlO'7  cm’3  (co-  736)° 764.  (14) 

In  the  above  approximation  <y  is  the  A\(LO)  phonon-plasmon  coupled  mode  in  cm'1,  and  the  value 
736  is  the  frequency  of  the  A\(LO)  Raman  mode.  At  that  low  level  of  free  carriers,  Ponce  et  al. 
investigated  the  spatial  distribution  in  GaN  crystallites  via  a  Raman  imaging  technique  [56].  The 
underlying  principle  of  the  technique  involves  the  recording  and  digitization  of  the 
/4i(TQ)-plasmon  mode  intensity  across  the  sample.  The  images  presented  by  Ponce  et  al.  show  a 
variation  of  brightness  across  the  crystallite,  with  the  areas  of  brightness  corresponding  to  low 
doping  concentration  [56]. 

7.  Isotopic  Effects  and  Phonon  Lifetimes  in  the  Wurtzite  Materials 

Since  the  isotope  mass  affects  the  phonon  frequency,  Raman  analysis  of  isotopic  films  may  convey 
information  on  the  elements  controlling  the  mode-vibrations.  Such  studies  of  WZ-GaN  films  made 
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from  natural  Ga  and  N  as  well  as  the  isotope  15N  have  been  reported  by  Zhang  et  al.  [63].  In  the 
film  containing  the  l5N  all  of  the  observed  Raman  modes,  excluding  the  E\ ,  were  found  to  exhibit 
prominent  downward  frequency  shifts.  Additionally,  the  Raman  frequencies  of  the  polar  modes 
Ax(LOJO)  and  E\(LO,TO)  were  found  by  Zhang  et  al.  to  be  shifted  according  to  the  inverse 
square  root  of  the  reduced  masses,  as  would  be  expected  from  the  first  principal  calculations  of 
crystal  dynamics  [63].  Since  the  reduced  mass  is  ll  —  l///7Ga  T  1  /v?n,  the  vibrations  of  the  polar 
modes  are  mainly  due  to  the  nitrogen  oscillations.  However,  the  isotope  shifts  of  the  non-polar 
modes  E\  and  E\  were  found  to  deviate  from  the  expected  reduced  mass  behavior,  with 
significant  deviation  found  for  the  E\  which  seems  not  to  be  affected  by  the  isotopic  mass.  The 
model  calculations  used  to  explain  the  vibrational  dynamics  of  the  non-polar  modes  led  the 
authors  to  conclude  that  although  the  E2  mode  is  still  dominated  by  nitrogen  atom  vibrations,  it 
constitutes  only  92%  of  the  total  vibrations.  The  rest  of  the  8%  vibrations  in  the  admixture  are 
due  to  the  movements  of  the  Ga  atoms.  A  different  scenario  was  inferred  for  the  E2  mode,  its 
vibration  mostly  involves  the  motion  of  Ga  which  explains  the  lack  of  frequency  response  to  the 
isotope  mass  of  the  ,5N.  Thus,  Zhang’s  studies  indicated  that  all  of  the  Raman  modes  but  the  E\ 
involve  the  vibrations  of  the  nitrogen  atoms,  while  that  of  the  E2  involves  the  vibrations  of  the 
heavier  gallium  atoms  [63]. 

One  crucial  aspect  impacting  device  performance  is  phonon  interaction  with  free  carriers.  In 
general,  the  interaction  can  degrade  the  viability  of  the  device;  however,  studies  have  also 
demonstrated  that  the  phonon-electron  interaction  may  be  used  to  engineer  certain  laser  devices 
[64-66].  The  phonon  lifetimes  are  important  in  both  these  aspects,  and  although  the  interaction 
involves  only  the  LO  phonons,  knowledge  of  other  mode-lifetimes  may  give  insight  into  the 
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characteristic  dynamics  of  the  material.  One  fundamental  lifetime  shortening  mechanism  in 
semiconductors  has  been  established  to  occur  via  the  anharmonic  interaction  [67-72].  In  this 
mechanism,  the  Raman  phonons  decay  into  other  normal  modes  in  such  a  way  that  there  is  a 
conservation  of  momentum  and  energy  in  the  process.  More  specifically,  for  a  three-phonon  decay 
process  a  phonon  of  frequency  <y,  and  a  wavevector  q\  decays  into  two  phonons  of  energies 
co 5,  and  wavevectors  <71,  qz,  such  that  d>  =  o>  ■+  cOj  and  q \  =  q_  +  (]}■  However,  the  additional 
lifetime  shortening  mechanism  due  to  phonon  scattering  at  point  defects  has  to  be  considered  as 
well  [67], 

Tsen  et  al.  have  reported  phonon  lifetime  measurements  via  Raman  spectroscopy  [73].  In  their 
investigation,  the  decay  of  the  Ai(LO)  mode  in  WZ-GaN  film  was  studied  via  time-resolved 
Raman  spectroscopy.  The  measured  lifetime  was  found  to  be  ~  3  ps  at  300  K  and  ~5  ps  at  5K. 
Tsen  et  al.  hypothesized  that  the  zone-center  A 1  (LO)  phonons  decay  primarily  into  a  large 
wavevector  TO  phonon  and  a  large  wavevector  LA  or  TA  phonon.  [73]. 

Raman  studies  of  phonon  lifetimes  in  GaN,  AIN,  and  ZnO  crystallites  have  been  reported  by 
Bergman  et  al.  [27].  The  lifetimes  were  obtained  from  the  Raman  linewidth,  after  correcting  for 

AE  1 

the  instrument  broadening,  using  the  uncertainty  relation  —  =  — ,  where  AE  is  the  linewidth  and 

tl  T 

ris  the  lifetime  [74].  The  lifetime  analysis  of  Bergman  et  al.  indicated  that  the  phonon  lifetimes  in 
AIN,  GaN,  as  well  as  ZnO  crystallites  fall  into  two  main  time  regimes:  a  relatively  long  lifetime  of 
the  .E\  mode  and  much  shorter  lifetimes  for  the  E\ ,  £i(TO),  A\(TO),  and  A\(LO)  modes.  The 
lifetime  of  the  E\  mode  of  high  quality  GaN  crystallites  was  found  to  be  «  10  ps,  whereas  the 
lifetimes  of  the  other  modes  were  found  to  be  approximately  an  order  of  magnitude  shorter.  A 
similar  trend  in  the  lifetimes  was  observed  for  phonons  of  high  quality  AIN,  ZnO,  and  AIN  that 
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contain  high  levels  of  impurities.  Moreover,  the  lifetimes  of  all  the  modes  of  both  the  high  and  low 
quality  AIN  crystallites  were  found  to  be  correlated  to  their  relative  impurity  concentrations  [27]. 

A  tentative  explanation  of  the  relative  long  lifetime  of  the  E\  in  the  WZ-crystallites  was  given 
by  Bergman  et  al.  in  terms  of  the  factors  determining  the  anharmonic  lifetimes,  specifically  the 
energy-conservation  constraints,  the  density  of  the  final  states,  as  well  as  the  anharmonic 
interaction  coefficient.  Unlike  the  other  modes,  the  energy  of  the  E\  mode  lies  in  the  low  energy 
regime  of  the  wurtzite  dispersion  curve  [75-77]  and  only  the  acoustic  phonons  are  available  as  a 
channel  of  decay.  At  the  zone  edges  the  energies  of  the  acoustic  phonons  are  equal  or  larger  than 
that  of  the  E\  mode;  thus,  in  order  for  the  energy  conservation  to  hold  the  optical  phonons  have 
to  decay  into  acoustic  phonons  at  the  zone-center  for  which  their  density  is  low.  Although  the 
contribution  of  the  anharmonic  coefficient,  which  has  not  yet  been  determined  theoretically,  has  to 
be  taken  into  account,  the  authors  have  suggested  that  the  low  density  of  states  significantly 
reduces  the  scattering  rate,  thus  increasing  the  phonon  lifetime. 

8.  Wide  Band-Gap  Alloys 

In  the  following,  a  review  of  the  mode  behavior  of  III-V  nitride-based  alloys  is  presented.  Mixed 
crystals  of  the  form  ABiX*  are  classified  into  two  main  groups  according  to  the  behavior  of  the  q 
«  0  optical  phonons  [78].  Figure  24  depicts  the  two  classes  for  ZB  material,  referred  to  as  the 
one-mode  and  two-mode  material,  respectively  [79].  In  general,  if  the  frequencies  of  the  AB  and 
the  AC  components  differ  greatly,  a  two-mode  behavior  is  expected;  if  the  frequencies  of  both 
components  have  proximate  values,  a  one-mode  behavior  results.  In  addition  to  the  one-mode  and 
two-mode  classes  of  materials,  an  intermediate  class  exists  that  exhibits  two-mode  behavior  over  a 
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The  polar  phonon  in  III-V  ternary  nitride  semiconductors  of  wurtzite  structure  were 
investigated  theoretically  by  Yu  et  al.  [86].  Within  the  modified-random-element  isodisplacement 
model  (MREI),  they  demonstrated  that  the  polar  modes  of  Ga^Al^N  and  ImGa,..rN  follow  a  one¬ 
mode  behavior  [86].  The  MREI  model  in  [86]  is  based  on  the  model  developed  previously  for  the 
ZB  materials  with  the  modifications  that  include  the  additional  phonon  modes  and  the  anisotropy 
of  the  WZ  structure.  Figures  25  and  26  present  the  theoretical  results  and  some  of  the  compiled 
experimental  results.  As  can  be  seen  in  Fig.  25,  the  theory  concurs  well  with  the  experiment 
except  for  the  Ex(TO)  mode.  This  small  deviation  has  been  attributed  to  crystal  quality.  The  one¬ 
mode  behavior  found  by  Yu  et  al.  was  explained  in  terms  of  the  large  mass  difference  between  the 
nitrogen  and  the  other  alloy  constituents.  Since  the  atomic  mass  of  nitrogen  is  much  smaller  that 
that  of  gallium,  indium,  and  aluminum,  the  reduced  masses  of  GaN,  InN,  and  AIN  are  almost  the 
same  as  that  of  nitrogen.  As  a  result,  no  distinct  GaN-like  modes  and  AlN-like  modes  may  exist  in 
the  Ga^AliJS!  alloy  system;  the  same  situation  holds  for  the  ImGa^N  alloy  [86],  It  should  be 
noted  that  in  all  of  the  aforementioned  experimental  results  the  other  possible  factors  affecting  the 
line  position  such  as  stress,  phonon-plasmon  coupling,  and  symmetry  mixing  were  not  taken  into 
consideration. 

Raman  studies  on  phonon  mode  behavior  in  ZB-AlGaN  films  were  reported  by  Harima  et  al. 
[15].  It  was  found  that  the  LO  mode  is  a  one-mode  type  while  the  TO  mode  is  a  two-mode  type. 
The  authors  calculated  the  mode  behavior  for  the  ZB-AlGaN  alloy  system  and  found  it  to  be 
consistent  with  the  experimental  results.  Figure  27  presents  the  phonon  frequency  versus  the 
composition  for  the  experimental  as  well  as  the  calculated  results  [15], 

The  topic  of  ordering  in  the  AlGaN  alloy  system  was  addressed  by  Korakakis  et  al.,  utilizing 
X-ray  diffraction  [87],  and  by  Bergman  et  al.,  utilizing  Raman  and  X-ray  spectroscopy  [88],  The 
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films  in  the  former  study  were  grown  by  MBE'at  a  temperature  of  «  750°C  on  SiC  and  sapphire 
substrates.  The  ordering  in  the  films  was  inferred  from  the  superlattice  (SL)  lines  of  the  (0001), 
(0003),  and  (0005)  diffraction  [87].  These  lines  are  forbidden  reflections  in  WZ  structure  and 
appear  as  SL  lines  only  when  lattice  ordering  exists.  The  intensity  ratio  of  the  SL  line  to  the 
allowed  (0002)  line  of  the  WZ  structure  is  thus  a  measure  of  the  extent  of  the  order.  The  X-ray 
analysis  of  Korakakis  indicated  a  long-range  order;  the  relative  intensity  of  SL  lines  was  found  to 
be  strong.  In  contrast  to  these  results  the  AlGaN  films,  which  were  grown  via  MOCVD  at  the 
elevated  temperature  range  »  1100°C  and  which  have  been  investigated  by  Bergman  et  al.,  were 
found  to  exhibit  a  much  smaller  intensity  ratio  thus  implying  a  disordered  state  of  the  alloy  [88]. 
The  different  order  states  found  in  both  studies  may  be  attributed  to  the  temperature  employed  in 
growth,  since  it  has  been  demonstrated  that  the  achievement  of  an  ordered  alloy  in  some  families 
of  ternary'  tetrahedral  semiconductors  is  a  function  of  the  growth  temperature  [89]. 

The  influence  of  the  disorder  on  the  E;  Raman  mode  of  the  ALGa^N  MOCVD  films  in  the 
composition  range  0<x<l  has  been  investigated  as  well  by  Bergman  et  al.  [88],  Figure  28  presents 
the  Raman  spectra  of  the  E\  mode  at  various  alloy  compositions;  it  is  shown  that  the  spectral 
lines  exhibit  asymmetric  broadening  toward  the  higher  frequency  range.  The  asymmetry  of  the  E\ 
lines  was  explained  in  light  of  the  spatial  correlation  model  [90-92].  The  foundation  of  the  model 
lies  in  the  wave  vector  uncertainty:  Aq=2n/L.  In  that  relation  Aq  is  the  phonon  wavevector  range; 
this  range  is  due  to  the  relaxation  of  the  q  =  0  Raman  selection  rules  caused  by  the  disorder.  The 
parameter  L  may  be  regarded  as  the  size  of  the  embedded  ordered  domain  in  the  disordered 
matrix:  thus  the  larger  the  disorder,  the  smaller  the  ordered  domain.  The  model  may  be  expressed 
as [90-92] 
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1{(0)  oc  jexp(-?2I2  /  4)— ^2  +  [fo  /  2j2  • 


(15) 


In  this  relation  77,  is  the  linewidth  of  the  alloy  of  composition  x  =  0  and  cc(q)  is  the  phonon 
dispersion  curve.  Few  and  to  some  extent  conflicting  theoretical  predictions  of  the  dispersion 
relation  in  the  III-V  nitrides  semiconductors  have  been  reported  [75,  76,  93-99].  Moreover,  due 
to  the  lack  of  experimental  results  no  consensus  has  been  reached  regarding  the  phonon  dispersion 
relations  in  these  materials.  In  order  to  fit  the  Raman  data  to  the  model,  the  authors  in  [88] 
assumed  an  averaged  value  of  co(q)  of  the  form  A+Bq2  (in  units  of  cm*1  and  with  ,4=568  and 
£=100).  As  can  be  seen  in  Fig.  28,  the  Raman  data  for  the  x  =  0.12  and  x  =  0.22  can  be  fit  with 
this  model  (represented  by  the  solid  lines),  albeit  with  different  values  for  77  These  fits  suggest  an 
estimate  I  «  3.5  nm  for  the  size  of  the  embedded  ordered  domain  [88].  Figure  29  shows  the  E2 
Raman  linewidth  as  a  function  of  composition;  a  maximum  is  observed  at  composition  x  »  0.5,  a 
value  at  which  maximum  disorder  would  be  expected  in  a  random  alloy. 

9.  Concluding  Remarks 

The  Raman  effect  arises  from  the  interaction  of  light  with  matter;  as  such,  Raman  spectroscopy  is 
a  nondestructive  and  powerful  tool  in  the  study  of  lattice  dynamics.  The  various  investigations 
reviewed  here  demonstrate  the  utility  of  Raman  spectroscopy  in  characterizing  the  material 
properties  of  the  III-V  nitrides.  Although  significant  progress  has  been  made  many  issues  are  still 
in  the  initial  stages  of  research.  The  cubic  phases,  the  InN  and  its  alloy  system,  and  the 
deconvolution  of  the  combined  phonon  dynamics  effect  on  the  alloy  Raman  spectra  merit  closer 
investigation.  Moreover,  the  topic  of  Raman  spectroscopy  of  phonon-interfaces  and  phonon- 
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superlattices  needs  to  be  addressed.  The  successful  investigation  of  many  of  the  above  issues  is 
contingent  on  advances  in  material  quality  of  the  wide  band-gap  semiconductors. 
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Figure  Captions 

Fig.  1 .  Optical  phonon  modes  in  the  wurtzite  structure.  The  Raman  active  modes  are  the  polar 
A  i  and  Ex ,  and  the  two  non-polar  E2  modes.  The  Bx  modes  are  silent  [19]. 

Fia.  2.  Schematic  phonon  dispersions  for  the  ZB  and  the  WZ  structures.  The  phonon  curves: 
TA,  LA,  TO,  and  LO  of  the  ZB  are  along  the  T-L  [111]  direction  ( thick  lines).  The 
folding  (thin  lines)  of  these  curves  into  the  WZ  Brillouin  zone  T-A  [0001]  direction 
creates  the:  E\  from  the  TA,  B\  from  the  LA,  Ej  from  TO  and  B{  from  LO.  The 
splitting  (dashed  lines)  creates  the  A\(TO)  and  the  Ex(TO)  from  the  original  TO  as  well 
as  the  Ax(LO)  and  the  Ex(LO)  from  the  original  LO. 

Fig.  3.  Polarized  Raman  spectra  of  GaN  film  in  the  v-direction  backscattering  geometry,  the 
GaN  modes  are  present  as  well  as  the  sapphire  substrate  modes  [26]. 

Fig.  4.  Polarized  Raman  spectra  of  GaN  for  the  z-direction  back- scattering  Raman  configuration 
[26]. 

Fig.  5.  Depolarized  Raman  spectra  of  GaN  films  grown  on  GaAs  substrates.  Sample  A  is  mainly 
WZ,  sample  C  is  mainly  ZB,  and  sample  B  is  a  mixed  phase  structure  [17]. 

Fig.  6.  Raman  spectra  of  WZ-GaN  film  at  three  incident  laser  wavelengths:  (a)  514.5  nm,  (b) 
488  nm,  and  (c)  457.9  nm  [14]. 

Fig.  7.  Raman  spectra  of  the  ZB-GaN  film  at  the  three  incident  laser  wavelengths  as  in  Fig.  6 

[14]. 

Fig.  8.  Raman  spectra  and  the  curve  fit  to  the  spectra  of  AIN  crystallite  [27]. 

Fig.  9.  Compositional  variation  of  the  TO-Raman  mode  of  ZB-AlGaN  alloy.  The  phonon 
frequency  of  the  ZB-AIN  is  655  cm  1  [15]. 

Fig.  10.  LO  phonon  spectra  of  ZB-AlGaN  ahoy.  The  LO  mode  of  the  ZB-AIN  is  observed  at  902 
cm'  [15]. 

Fig.  1 1 .  Model  calculations  (lines)  of  the  high  frequency  phonon  of  WZ-AIN  and  GaN  under 
hydrostatic  pressure.  The  dots  represent  experimental  data  from  [30,  38, 19]. 

Fig.  12.  Model  calculations  (lines)  of  the  low  frequency  phonons  of  WZ-AIN  and  GaN  under 
hydrostatic  pressure.  The  dots  represent  experimental  data  from  [30,  38, 19]. 

Fig.  13.  Raman  frequency  scheme  for  (a)  the  case  when  the  long  range  electrostatic  field  is  the 
dominant  interaction,  and  (b)  when  the  anisotropy  of  the  short  range  interaction 
dominates  [49]. 
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Fig.  14.  Geometry  of  the  Raman  experiment  setup  for  the  observation  of  the  quasi-modes.  The 
spectra  were  acquired  in  a  backscattering  geometry  from  the  Y  axis,  6  is  the  rotation 
angle,  and  p  is  the  angle  of  the  phonon  propagation  [49]. 

Fig.  15.  Raman  spectra  at  various  rotations  and  Raman  configuration  which  enable  the 
observation  of  the  quasi- TO  modes.  The  pure-modes:  E] ,  and  E{(TO)  are  present  as 
well  [49]. 

Fig.  1 6.  Raman  spectra  of  the  quasi-LO  modes  of  an  AIN  crystallite  [49]. 

Fia.  17.  Angular  dispersion  of  the  quasi-IO  and  -TO  modes.  The  lines  represent  Loudon’s 
model;  the  dots  and  Xns  are  the  experimental  results  from  [49]  and  [32],  respectively 
[49],  ’ 

Fig.  18.  Raman  spectra  of  the  A\(LO)  mode  of  GaN  at  several  Si  doping  concentrations  nH 
obtained  via  Hall  measurements  [51]. 

Fig.  19.  Experimental  and  calculated  Raman  lineshape  of  the  ^iflQ)-plasmon  coupled  mode 

[51]- 

Fig.  20.  Correlation  between  the  carrier  concentration  n  found  from  the  Raman  analysis  and  the 
concentration  n#  obtained  via  Hall  measurements  [51]. 

Fig.  21.  Raman  spectra  of  highly  conductive  GaN  films.  Dotted  line :  spectra  from  sample  with 
free  carrier  density  n  =  4xl019  cm'3;  solid  line:  from  sample  with  n  =  8xl019  cm'3  [54]. 

Fig.  22.  Calculated  spectra  of  the  coupled  plasmon-phonon  modes  for  a  sample  with  n  =  4xl019 
cm'3.  The  dotted  line  corresponds  to  a  plasma  damping  constant  of  100  cm' ;  the  solid 
line  to  1 1 16  cm'1,  which  is  a  realistic  value  for  the  GaN  film  [54], 

Fig.  23.  Experimental  and  calculated  ( dotted  line )  lineshape  of  the  low-branch  A\(LO)  phonon- 
plasmon  coupled  mode  in  GaN  films  [55]. 

Fig.  24.  The  two  classes  of  mixed  ZB  crystals:  (a)  two-mode  type  and  (b)  one-mode  type 
material  [79]. 

Fig.  25 .  Phonon  mode  behavior  of  WZ-GaAIN  alloy.  The  lines  represent  the  MRE1  model.  The 
model  calculations  are  for  the  A\(TO),  A\(LO),  E\(TO),  and  E\(LO)  modes;  the  points 
represent  the  experimental  data  [86]. 

Fig.  26.  Phonon  mode  behavior  of  WZ-InGaN  alloy.  Theory  (lines')  and  experiment  (dots).  The 
polar  modes  are  presented  [86]. 

Fig.  27.  Mode  behavior  of  ZB-AlGaN.  For  an  intermediate  composition  x  in  the  two-mode  class 
of  materials,  two  sets  of  frequencies  would  be  observed  in  the  spectra:  one  set  is  due  to 
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the  LO  and  TO  modes  of  the  AB  component  and  the  other  set  is  due  to  the  LO  and  TO 
modes  of  the  AC  component.  The  one-mode  exhibits  only  one  set  of  LO  and  TO 
frequencies,  which  ideally  are  linear  with  composition  [15]. 

Fig.  28.  Raman  spectra  of  the  E\  mode  of  WZ-AlGaN  films  at  various  compositions.  The  solid 
lines  represent  the  spatial  correlation  model  [88]. 

Fig.  29.  The  E]  Raman  linewidth  as  a  function  of  composition  (dots).  The  solid  line  is  obtained 
from  calculations  of  the  entropy  of  mixing  of  an  alloy  system  (in  arbitrary  units)  [88]. 
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We  report  on  UV  Raman  spectroscopy  of  InxGai.xN  thin  films  grown  on  sapphire 
by  metal-organic  chemical  vapor  deposition  (MOCVD).  For  the  first  time  the  Ai(LO) 
and  E2  phonon  mode  behavior  was  investigated  over  a  large  compositional  range 
(0<x<0.50).  Compelling  evidence  is  presented  for  one-mode  behavior  for  the  A](LO) 
phonon,  and  data  suggestive  of  two-mode  behavior  is  presented  for  the  E2  phonon. 
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In  recent  years  interest  has  increased  in  the  group  III  nitrides  as  hopeful 
candidates  for  the  fabrication  of  light  emitting  diodes  and  laser  diodes  operating 
throughout  the  visible  region  and  into  the  UV  range.  AIN,  GaN,  and  InN  are  an 
important  set  of  semiconductors  because  their  direct  band  gaps  span  the  range  1 .9  to 
6.2eV  [1].  In  principle,  this  complete  set  of  ternary  alloys  allows  for  any  band  gap  within 
this  range.  This  tunability  of  the  energy  band  gap,  combined  with  the  fact  that  they 
generate  efficient  luminescence,  has  motivated  the  study  and  fabrication  of  high 
performance  optical  devices  based  on  the  group  III  nitrides.  Potential  applications  range 
from  high-density  optical  storage  devices  to  projection  and  other  bright  displays. 
Presently  much  work  and  progress  has  been  done  on  AIN,  GaN,  and  their  alloys.  The 
growth  of  InGaN  alloys,  with  a  band  gap  energy  range  of  1.9  to  3.4eV,  however,  has 
proven  to  be  very  difficult.  One  of  the  main  difficulties  with  InGaN  growth  is  the  large 
difference  in  interatomic  spacing  between  GaN  and  InN  of  roughly  1 1%.  This  large 
difference  can  give  rise  to  a  solid  miscibility  gap  and  highly  strained  InGaN  alloys  [2], 

As  a  result,  phase  separation  in  the  alloy  becomes  a  major  concern.  Until  recently,  there 
has  been  little  attention  given  to  the  phonon  mode  behavior  over  a  wide  compositional 
range  in  the  hexagonal  wurzite  InxGai.xN  system. 

Mid  and  deep-UV  Raman  spectroscopy  has  been  used  here,  we  believe  for  the 
first  time,  as  an  effective  tool  for  investigating  LO-plasmon  interactions  and  studying  the 
issue  of  compositional  inhomogeneities  and  spinodal  decomposition  in  InGaN  thin  films. 
In  addition,  mode  behavior  over  a  large  compositional  range  was  investigated,  and  we 
have  found  compelling  evidence  for  one-mode  behavior  for  the  Ai(LO)  phonon  mode. 


while  our  data  is  suggestive  of  two-mode  behavior  for  the  E2  mode  consistent  with  the 
previously  found  behavior  in  AlGaN  [3,4]. 

In  this  work,  we  study  the  Raman  scattering  spectra  of  InxGai.xN  alloys,  grown  on 
(0001)  sapphire  substrates  using  a  specially  designed  metal-organic  chemical  vapor 
deposition  (MOCVD)  reactor  and  growth  conditions  that  have  been  described  elsewhere 
[5],  A  buffer  layer  of  GaN  was  grown  on  each  substrate  prior  to  the  growth  of  the  InGaN 
films.  The  growth  temperature  of  the  InxGai.xN  layer  was  between  690°C  and  780°C. 
Variation  of  the  growth  temperature  and  the  hydrogen  gas  flow  rate  determined  the 
percent  concentration  of  InN.  Eight  films  were  examined  in  the  compositional  range 
0<x<0.50.  The  composition  x  was  determined  by  calculating  the  shift  of  the  0002  XRD 
InGaN  peak  relative  to  the  0002  peak  position  of  GaN,  and  applying  Vegard’s  law. 

Raman  scattering  experiments,  of  thin  InxGai_xN  films,  using  visible  laser 
frequencies  typically  provide  additional  Raman  peaks  from  buffer  layers  and  substrates. 

In  addition,  InGaN  films  with  low  InN  compositional  content,  emit  strong  interfering 
broadband  luminescence  in  visible  regions.  The  excitation  of  Raman  scattering  in  the 
UV  avoids  both  of  these  problems.  The  Raman  spectral  region  of  interest  is  well  above 
any  emitted  photoluminescence,  and  there  is  strong  absorption  of  the  UV  laser  light  in  the 
film,  resulting  in  a  larger  Raman  signal  and  an  absence  of  Raman  peaks  from  layers 
below  the  film.  The  deep-UV  Raman  spectra  were  excited  with  the  244nm  line  of  a 
frequency  doubled  Coherent  Inc.  Innova  300  FreD  Argon  Ion  laser.  The  incident  beam 
was  focused  onto  the  sample  to  a  spot  size  of  approximately  2mm  in  height  and  50 
microns  in  width.  The  intensity  of  the  light  on  the  surface  of  the  sample,  in  this 
configuration,  was  on  the  order  of  65  W/cm2.  The  scattered  light  was  collected  in  a 


backscattering  geometry  and  imaged  onto  the  entrance  slit  of  a  SPEX  1877  TripleMate 
spectrometer,  equipped  with  a  UV  enhanced  liquid  nitrogen  cooled  CCD.  Specialized 
high  reflectance  dielectric  mirrors  were  used  for  steering  the  laser  beam,  and  a  single 
piano  convex  lens,  of  UV  grade  fused  silica,  was  used  to  focus  the  scattered  light.  The 
Raman  spectra  were  obtained  with  a  3600  gr/mm  grating,  providing  0.8  cm'1  resolution. 
The  mid-UV  Raman  spectra  were  excited  with  the  325  nm  line  of  a  HeCd  laser.  The 
laser  light  was  focused  on  the  sample  through  a  specialized  lOOx  UV  microscope 
objective  down  to  a  spot  size  of  approximately  0.3  microns  with  a  power  density  around 
1 6.7x1 05  W/cm2.  The  scattered  light  was  analyzed  with  a  UV  JY  LabRam  Infinity 
equipped  with  a  UV  enhanced  liquid  nitrogen  cooled  CCD.  The  Raman  spectra  were 
obtained  with  a  2400  gr/mm  grating  providing  3  cm'1  resolution.  For  both  of  these 
situations  Raman  scattering  by  the  Ai(LO)  and  E2  phonons  are  symmetry  allowed  from 
the  (0001)  surfaces  of  the  hexagonal  wurzite  structures. 

Figure  1  shows  Raman  spectra  excited  with  the  244nm  (5.1  eV)  line  of  the 
frequency  doubled  Ar  ion  laser.  The  dashed  lines  indicate  the  Raman  frequencies  of  the 
E2  mode  for  InN  and  GaN.  The  value  for  the  E2  mode  in  InN  was  taken  from  other 
references  [3-6].  Several  Raman  peaks  due  to  scattering  from  air  overlap  the  region  of 
the  Ai(LO)  peak.  Analysis  of  the  A](LO)  peak  positions,  for  the  data  taken  with  the 
244nm  line,  was  not  included  for  this  reason.  However,  qualitative  observations 
regarding  the  ratio  of  Ai(LO)  to  E2  peak  intensities  are  made.  For  the  GaN  and  the  films 
with  low  InN  concentration,  the  intensity  of  the  Ai(LO)  mode,  with  respect  to  the  E2 
mode,  is  small.  The  absence  of  the  Ai(LO)  mode  has  been  previously  observed  in  GaN 
and  has  been  attributed  to  a  high  free  carrier  concentration  in  the  film;  the  free  carriers 
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form  plasmon  waves  that  interact  and,  in  turn,  damp  the  LO  phonons  [10,1 1],  The  free 
carrier  concentration  of  the  GaN  film  used  for  Figure  1  was  found  to  be  8.4xl016  cm'3  by 
Hall  measurements.  The  free  carrier  concentrations  in  the  InGaN  films  are  of  the  order  of 
1018  cm'1,  and  the  concentration  increases  with  higher  In  concentration.  Despite  this,  the 
films  with  higher  InN  concentration  show  an  increase  in  the  Ai(LO)  intensity  with 
respect  to  the  E2  mode.  Inushima  et  al.  reported  Raman  spectra  of  InN  films  in  which  the 
Ai(LO)  mode  intensity  was  also  greater  than  the  E2  mode  intensity  [9].  Their  films  were 
grown  on  (0001)  sapphire  substrates  with  a  carrier  density  of  3x1 020  cm"J  by  Hall 
measurements.  They  ascribed  this  effect  to  the  presence  of  a  deep  depletion  layer  on  the 
air  exposed  surface  of  InN. 

Figure  2  shows  Raman  spectra  excited  with  the  325. 2nm  (3.8  leV)  line  of  the 
HeCd  laser.  The  dashed  lines  again  show  the  position  of  the  E2  mode  for  InN  and  GaN. 
Raman  spectra  for  GaN  could  not  be  obtained  due  to  the  strong  photoluminescence 
resulting  from  excitation  just  above  the  band  gap  of  the  material.  In  these  spectra  the 
opposite  case  is  observed  than  in  Figure  1 .  The  E2  mode  intensity  is  weak  with  respect  to 
the  Ai(LO)  mode  for  low  concentrations  of  InN  in  the  film.  The  Ai(LO)  mode  is  easily 
detected,  and  though  not  shown  here,  even  second  and  third  order  Ai(LO)  phonons  could 
be  observed.  This  is  presumably  due  to  a  resonant  Raman  effect  at  this  laser  energy. 

Also  worth  noting,  the  ARLO)  mode  develops  an  additional  low  energy  feature  for  In 
concentrations  higher  than  15%.  The  shoulder  and  the  large  FWHM  of  the  ARLO)  peaks 
in  Figure  2  have  tentatively  been  assigned  to  compositional  inhomogeneities  and  spinodal 
decomposition  effects  present  in  InGaN  films  that  have  been  previously  discussed 
[2,12,13].  This  issue  will  be  addressed  more  completely  in  a  future  publication. 
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Figure  3  shows  the  frequencies  of  the  Ai(LO)  and  E2  modes  as  a  function  of  alloy 
composition.  The  filled  circles  represent  the  peak  positions  of  the  Ai(LO)  mode  in 
Figure  2.  The  peak  positions  of  the  Ai(LO)  mode,  in  Figure  1,  are  not  shown  due  to  the 
presence  of  interfering  air  peaks,  which  made  assignment  of  peak  positions  difficult.  The 
crosses  and  unfilled  circles  in  Figure  3  represent  the  E2  mode  frequencies  from  Figures  1 
and  2  respectively.  The  peak  positions,  for  the  data  in  Figures  1  and  2,  were  determined 
by  fitting  the  spectra  with  Lorentzian-Gaussian  peak  functions.  The  Ai(LO)  peaks  in 
Figure  2,  which  exhibited  a  low  energy  shoulder,  were  fit  with  two  Lorentzian-Gaussian 
peak  functions.  The  sample  with  35.6%  In  concentration  was  fit  with  an  additional  high 
energy  shoulder.  The  Ai(LO)  peak  positions  used  in  Figure  3  are  marked  in  Figure  2. 
These  positions  are  assumed  to  be  due  to  the  majority  phase  of  the  alloy  in  the  films,  and 
the  low  energy  shoulder  is  tentatively  assigned  to  the  phase  segregated  regions  of  higher 
In  concentration.  Similar  InGaN  films  were  analyzed  by  0-20  x-ray  diffraction.  The  x- 
ray  data  for  films  with  concentrations  above  40%  In  concentration  showed  a  diffraction 
peak  for  the  majority  phase  along  with  a  shoulder  that  was  fitted  with  several  Gaussian 
distribution  functions  [13].  The  authors  ascribed  this  feature  to  phase  separation 
occurring  in  the  films.  The  value  for  the  Ai(LO)  mode  in  InN  was  taken  from  other 
references  [6-9,14],  The  Ai(LO)  frequencies  shift  linearly  as  a  function  of  percent  InN 
composition.  The  experimental  values  presented  here  fit  well  with  a  linear  extrapolation 
to  the  InN  endpoint.  This  is  indicative  of  one-mode  behavior  for  the  Ai(LO)  mode.  The 
one-mode  behavior  of  the  Ai(LO)  phonon  is  in  good  agreement  with  theoretical 
calculations  based  on  a  modified  random  element  isodisplacement  model  (MREI)  for  the 
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ternary  family  of  semiconductors  AB].XCX  [15].  It  has  also  been  shown  experimentally 
that  the  Ai(LO)  phonon  in  AlxGai.xN  wurzite  structures  obeys  a  one-mode  behavior  [3,4]. 

The  E2  phonon  frequencies  shown  in  Figure  3,  for  the  244nm  and  325.2nm  laser 
lines,  show  a  good  consistent  linear  fit  for  the  InN  compositional  range  of  0<x<0.50. 
However,  the  extrapolated  line  for  the  E2  data  intersects  between  30  to  40  cm'1  lower  at 
x=l  than  the  previously  reported  E2  mode  frequencies  for  InN.  It  has  been  shown 
experimentally  that  the  E2  phonon  mode  in  the  AlxGai_xN  wurzite  structure  displays  a 
two-mode  type  behavior  [3.4],  The  phonon  behavior  in  a  mixed  crystal  of  an  ABi.xCx 
ternary  system  falls  into  two  main  classes,  referred  to  as  one  mode  and  two  mode 
behavior.  Phonons  exhibiting  one  mode  behavior  have  only  one  set  of  frequencies  that 
show  nearly  linear  dependency  with  the  composition  x,  from  the  AB  component  to  the 
AC  component.  The  Ai(LO)  mode  in  InxGai.xN  exhibits  this  behavior.  On  the  other 
hand,  the  frequencies,  of  phonons  in  the  two  mode  class  of  materials,  exhibit  a  nonlinear 
dependence  with  composition  x.  Also,  for  an  intermediate  composition  x  in  the  two 
mode  class  of  materials,  two  sets  of  frequencies  are  usually  observed  in  the  spectra.  One 
set  of  frequencies  is  due  to  the  modes  of  the  lighter  AB  component  and  one  set  for  the 
heavier  AC  component.  These  components  in  the  InxGai.xN  crystal  are  GaN  and  InN. 

The  absence  of  the  InN  set  of  frequencies  in  the  spectra  may  be  due  to  the  large  FWHM 
for  peaks  close  to  50%  In  concentration,  making  it  difficult  to  detect  weak  peaks  in  that 
spectral  region. 

In  summary,  we  present  the  results  of  a  Raman  scattering  study  for  a  series  of 
InxGai_xN  films  from  x=0  to  0.50.  where  the  Raman  scattering  was  excited  with  two 
different  UV  frequencies.  Our  results  show  the  behavior  of  the  A|(LO)  and  E2  modes  as 
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a  function  of  frequency.  The  Ai(LO)  mode  follows  a  one  mode  behavior  while  the  E2 
mode  appears  to  follow  two  mode  behavior  consistent  with  AlGaN  studies.  The  two 
different  excitation  frequencies  provide  complimentary  information  as  to  the  aspects  of 
the  films,  and  provide  more  information  to  the  behavior  of  the  phonon-plasmon 
interaction  as  a  function  of  Indium  concentration  that  will  be  studied  further. 
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Figure  Captions 

Figure  1 :  Raman  Spectra  of  InGaN  films.  Spectra  excited  with 

the  244nm  line  of  the  frequency  doubled  Ar  ion  laser. 
The  dotted  lines  represent  the  positions  of  the  E2  mode 
for  InN  and  GaN  at  488  cm'1  and  738  cm'1. 

Figure  2:  Raman  Spectra  of  InGaN  films.  Spectra  excited  with 

the  325. 2nm  line  of  a  HeCd  laser.  The  dotted  lines 
represent  the  positions  of  the  E2  mode  for  InN 
and  GaN  at  488  cm’1  and  738cm'1.  The  positions  of  the 
Ai(LO)  frequencies  are  marked. 

Figure  3:  Percent  InN  concentration  vs  frequency  for  the 

Ai(LO)  and  E2  modes.  Ai(LO)  mode  for  325.2nm 
data.  E2  mode  for  both  244nm  and  325.2nm.  The  bars 
at  x=l  represent  data  points  for  InN  taken  from  the 
following  refs.  In  descending  order  [6,8,7,9,14]  for  the 
Ai(LO)  mode,  and  [6, 8, 7, 9]  for  the  E2  mode. 
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ABSTRACT 

Raman  analyses  of  the  lifetimes  of  phonons  in  GaN  and  AIN  crystallites  of  wurtzite  structure  are 
presented.  In  order  to  ensure  the  accuracy  of  the  measurement  of  the  phonon  lifetimes,  an 
experimental  procedure  to  eliminate  the  broadening  due  to  the  finite  slit  width  was  performed. 
The  lifetime  analyses  indicate  that  the  phonon  lifetimes  in  AIN  as  well  as  in  GaN  fall  into  two 
main  time  regimes:  a  relatively  long  time  of  the  E^  mode  and  much  shorter  times  of  the  E2 , 

El(TO),  and  Al(TO)  modes.  The  lifetimes  of  the  E^,  E;,  El(TO),  Al(TO),  and  Al(LO) 

modes  of  a  high  quality'  AIN  crystallite  are  4.4  ps,  0.83  ps.  0.91  ps,  0.76  ps,  and  0.45  ps, 
respectively.  Moreover,  the  lifetime  of  the  Al(LO)  mode  found  in  this  study  is  consistent  with 
the  current  phonon-decay  model  of  that  mode  in  wurtzite  structure  materials.  The  lifetimes  of 
E^ ,  E2 ,  El  (TO),  and  Al(TO)  of  a  GaN  crystallite  were  found  to  be  10.1  ps,  1.4  ps,  0.95  ps,  and 
0.46  ps,  respectively.  The  Al(LO)  mode  in  the  GaN  was  not  observed  and  its  absence  is 
attributed  to  the  plasmon  damping.  The  lifetime  shortening  due  to  impurities  was  also  studied; 
the  lifetimes  of  the  Raman  modes  of  an  AIN  crystallite  which  contains  about  two  orders  of 
magnitude  more  Si  and  C  impurities  relative  to  the  concentration  of  the  high  quality  crystallite 
were  found  to  be  50  %  shorter. 

PACS  Nos.:  63.20.-e,  63.20.Kr,  78.30.-j,  78.30.Fs 
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I.  INTRODUCTION 

The  realization  of  device  fabrication  using  AIN,  GaN,  and  their  alloy  systems  as  the 
constituent  materials,  has  led  to  numerous  investigations  concerning  growth  as  well  as  analytical 
techniques  to  study  and  improve  these  wide-bandgap  family  semiconductors  [1].  Raman 
spectroscopy  has  proven  useful  in  analyzing  the  properties  of  wide-bandgap  materials,  it  reveals 
information  about  sample  quality  as  well  as  phonon  interactions  and  dynamics  [2-10]. 

One  crucial  aspect  impacting  device  performance  is  phonon  interaction  with  free  carriers. 
In  general  the  interactions  can  degrade  the  viability  of  the  device;  however,  recent  studies  have 
demonstrated  that  the  phonon  interactions  may  be  utilized  to  engineer  certain  laser  devices  [11- 
12].  The  phonon  lifetimes  are  important  in  both  these  effects,  and  one  principal  way  of 
measuring  these  lifetimes  is  via  Raman  spectroscopy.  It  has  been  demonstrated  previously  that 
phonon  lifetimes  and  the  factors  affecting  them  can  be  calculated  from  the  Raman  linewidth. 
Among  the  materials  which  have  been  analyzed  via  their  Raman  linewidth  are  diamond,  Si,  c- 
BN,  GaAs,  and  InP  [13-17]. 

In  general,  the  theory  of  a  spectral  lineshape  of  a  signal  in  a  dispersive  medium  predicts 
the  lineshape  to  be  Lorentzian  and  the  linewidth,  a  parameter  describing  the  damping  effect,  to  be 
inversely  proportional  to  the  lifetime  of  the  signal  [18],  For  the  case  of  a  Raman  signal  of  an 
ideal  harmonic  crystal,  the  lineshape  is  expected  to  be  infinitesimally  narrow;  however, 
experiments  have  demonstrated  that  the  Raman  linewidths  of  most  materials  exhibit  a  finite 
width  indicative  of  the  presence  of  decay  channels  which  shorten  the  phonons  lifetimes  [19-21]. 
One  fundamental  lifetime  shortening  mechanism  in  semiconductors  with  of  diamond  structure 
has  been  established  to  occur  via  the  anharmonic  interaction  [20].  In  this  mechanism,  the  Raman 
phonons  decay  into  other  phonons  that  conserve  momentum  and  energy.  Moreover,  in 
semiconductors  containing  impurities  and  defects  it  has  been  found  that  the  these  imperfections 
affect  the  Raman  linewidth;  thus,  the  contribution  of  the  impurity  shortening  lifetime  mechanism 
has  to  be  taken  into  account  [21-22]. 

In  this  paper,  we  present  the  study  of  the  room-temperature  lifetimes  of  the  Raman 
phonons  in  AIN  and  GaN  materials.  The  Raman  modes  investigated  in  the  study  are  the  E^ , 
g 2  ?  El(TO),  Al(TO),  and  Al(LO)  modes.  The  Al(LO)  mode  of  the  GaN  crystallite  was  not 
observed  in  the  Raman  spectra  and  its  absence  is  attributed  to  the  high  free  earner  concentration 
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of  that  crystallite;  the  free  carriers  form  plasmon-waves  which  interact  and  in  turn  damp  the  LO 
phonons  [2,  7].  However,  due  to  the  insulating  properties  of  the  AIN.  the  Al(LO)  mode  was 
observed  in  the  spectra  of  the  crystallites.  For  AIN  and  GaN  crystallites  we  found  that  the 
phonons  can  be  separated  into  two  time  regimes:  a  long  lifetime  range  of  the  e\  phonon,  and  a 
much  shorter  lifetime  range  which  characterizez  phonons  including  the  E*,  El  (TO),  and  the 
Al(TO)  modes.  Moreover,  the  Al(LO)  phonon  of  the  AIN  was  found  to  have  the  shortest  value 
among  the  AIN  modes.  In  addition,  we  investigated  two  AIN  crystallites  which  differ  in  their 
impurity-concentration  and  found  that  the  phonon  lifetimes  were  correlated  with  the  impurity 
concentrations.  Also,  we  studied  the  E;  phonon  lifetime  of  high  quality  GaN  film  at  room 
temperature  and  at  10  K.  Our  results  indicate  that  the  phonon  lifetime  of  the  film  is  similar  to 
that  of  the  crystallite;  moreover,  no  significant  difference  between  the  lifetimes  in  the  two 
temperature  ranges  has  been  observed. 

The  Al(LO)  mode  lifetime  has  been  considered  theoretically  for  wurtzite-based  structures 
[23];  this  model  predicts  a  three-phonon  anharmonic  process  in  which  the  Al(LO)  decays  into  an 
El  (TO)  phonon  and  an  LA  phonon.  Our  experimental  measurements  of  this  lifetime  are 
consistent  with  the  predictions  of  the  model. 

II.  EXPERIMENT 

Single  crystals  of  GaN  and  AIN  were  grown  via  the  sublimation  method  described  in 
detail  in  [24-26].  A  short  summary  is  outlined  herein  highlights  the  findings  presented  in  these 
references  which  are  relevant  to  our  present  studies.  The  GaN  hexagonal  crystallites  are 
colorless  and  secondary  ion  mass  spectroscopy  (SIMS)  indicated  that  the  concentration  of  all 
impurities,  with  the  exception  of  oxygen,  were  at  background  levels;  the  oxygen  concentration 
was  3xl018  atoms/cm3.  Moreover,  photoluminescence  and  transmission  spectroscopies  indicated 
that  the  GaN  crystallites  were  of  high  optical  quality  and  no  deep  level  impurities  were  present. 

Two  different  qualities  of  hexagonal  AIN  crystallites  were  also  grown  via  the  sublimation 
method  [24]:  one  type  of  crystallites  was  grown  in  the  2100-2250  °C  temperature  range,  of  blue 
color,  and  will  be  referred  to  as  such  in  this  paper.  The  other  type  of  crystallites,  that  will  be 
referred  as  the  transparent  crystallites,  were  grown  in  the  1950-2050  °C  range.  SIMS  analysis 
confirmed  that  the  transparent  crystallites  contain  about  two  orders  of  magnitude  less  Si  and  C 
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impurities  than  the  blue  AIN.  Oxygen  impurity  concentrations  are  similar  in  both  samples.  The 
Si  and  C  impurities  were  incorporated  into  the  crystals  from  the  SiC  seed  substrates  while  the 
oxygen  was  incorporated  from  the  growth-environment.  The  transmission  electron  microscopy 
(TEM)  images  of  the  AIN  crystallites  indicated  that  they  are  of  high  quality;  no  high  angle 
boundaries,  stacking  faults,  or  twinned  region  were  observed.  In  addition.  X-ray  diffraction 
indicated  that  the  residual  stress  level  in  the  crystallites  is  low.  Fig.  1  depicts  one  of  the 
transparent  AIN  crystallites. 

The  micro-Raman  scattering  experiments  were  earned  out  at  room  temperature  utilizing 
the  514.5  nm  line  of  the  argon  ion  continuous  laser,  and  the  J-Y  U1000  scanning  double 
monochromator.  The  spot  size  of  the  laser  on  the  samples  was  ~  2  pm  in  diameter.  The 
experimental  error  of  the  Raman  data  is  ±  0.2  cm'1.  The  data  were  acquired  in  backscattering 

geometry  from  the  a-face  of  the  AIN  and  GaN  crystallites  for  the  observation  the  e\,  E* , 
El  (TO),  and  Al(TO)  modes,  and  from  the  c-face  for  the  observation  of  the  Al(LO)  mode.  The 
cold  temperature  measurements  were  earned  out  at  the  macro-mode  with  a  spot  size  ~T  mm  in 
diameter.  The  error  in  the  lifetime  values,  which  is  partially  due  to  the  expenmental  error  and 
partially  due  to  the  lineshape  data  fit,  was  estimated  to  be  ±  10%. 

HI.  RESULTS  AND  DISCUSSION 

In  order  to  obtain  the  phonon  lifetime  via  the  Raman  spectral  linewidth,  the  linewidth  has 
to  be  corrected  for  the  contribution  of  the  instrumental  bandpass  broadening.  The  following 
section  describes  the  method  used  in  this  study  to  determine  the  actual  phonon  linewidth.  The 
Raman  spectra  acquired  at  successive  slit  widths  ranging  from  400  pm  down  to  100  pm;  the 
measured  linewidths.  WM,  were  plotted  as  functions  of  the  slit  width,  Ws,  and  extrapolated  to  the 
zero-slit  value,  Wp,  via  the  relation  [27-28]: 

WM"iWP+(WS*9-2*10"3)  J» 

The  second  term  in  Equation  1  is  the  instrumental  bandpass,  i.e.,  the  slit  width  multiplied  by  the 
monochromator  linear  dispersion,  9.2*1 0‘3  ern'/pm.  This  method  yields  the  actual  phonon 
linewidth.  WP,  from  which  the  phonon  lifetime  may  be  calculated.  However,  Equation  1  is  an 
approximation  that  can  be  used  only  when  the  lineshape  of  the  instrument  bandpass  can  be 
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approximated  by  a  Gaussian  and  that  of  the  phonon  by  a  Lorentzian  [27-28],  To  check  the 
bandpass  characteristics  of  our  Raman  system,  a  spectrum  of  a  high  quality  diamond  crystal  (first 
checked  to  be  mostly  Lorentzian  at  100  pm  slits  and  Wp~  1.8  cm'1)  was  acquired  at  a  relatively 
open  slit  width  of  600  pm  in  which  case  the  dominance  of  the  instrumental  broadening 
component  is  expected.  The  lineshape  was  found  to  be  mostly  Gaussian  of  linewidth  5.5  cm 
which  can  be  attributed  to  instrumental  broadening.  Moreover,  the  lineshapes  of  the  nitride 
crystals  studied  here  exhibit  mainly  a  Lorentzian  component  at  the  lower  slit-values  as  depicted 
in  Fig.  2.  These  experimental  results  validate  the  appropriateness  of  using  Equation  1. 

In  general,  the  zero-slit  Raman  linewidth  arises  from  a  convolution  of  all  broadening 
mechanisms  which  in  turn  are  due  to  the  lifetime  shortening  mechanisms  occurring  in  a  given 
crystal  [21].  The  two  main  mechanisms  controlling  the  phonon  lifetimes  are  the  phonon 
anharmonic  interactions  in  which  a  phonon  decays  into  other  phonons,  and  the  phonon  scattering 
at  impurity  or  defect  centers  which  is  sometimes  treated  as  a  source  of  the  inhomogeneous 
broadening.  In  the  subsequent  discussion,  we  present  a  study  of  lifetime  shortening  mechanisms 
in  two  types  of  AIN  crystallites  which  differ  mainly  in  their  impurity  contents.  Figures  3  and  4 
depict  the  Raman  spectra  along  with  the  deconvoluted  spectra  of  the  AIN  crystallites.  The 
spectra  presented  in  Fig.  3  were  acquired  from  a  transparent  AIN  crystallite  while  those  in  Fig.  4 
are  from  a  blue  AIN  crystallite.  As  was  indicated  in  the  previous  discussion  on  experimental 
measurements,  the  blue  crystallite  contains  about  two  orders  of  magnitude  more  Si  and  C 
impurities  then  the  transparent  AIN.  In  both  spectra  the  Al(TO),  E; ,  and  El(TO)  Raman  modes 
are  at  608  cm'1,  655  cm'1,  and  668  cm'1,  respectively.  No  peak  shift  difference  between  the  two 
spectra  which  exceeds  the  experimental  error  of  ±  0.2  cm'1  is  evident.  Thus,  both  samples  are 
similar  in  their  high  quality  of  crystallinity  and  low  levels  of  residual  stress  state.  Moreover,  the 

impurities  have  no  significant  effect  on  the  residual  stress.  Figure  5  depicts  the  El2  Raman  mode 

at  246  cm'1  for  both  crystallites,  and  Fig.  6  presents  the  spectra  of  the  Al(LO)  mode  at  890  cm  . 

In  order  to  calculate  the  phonons  lifetime  via  the  Raman  linewidths,  the  linewidths  were 
plotted  as  a  function  of  the  slit  width  and  a  curve  fit  was  obtained  via  Equation  1.  The  results  are 
presented  in  Fig.  7;  the  dots  represent  the  data  taken  from  the  transparent  AIN  and  the  squares 
data  from  the  blue  crystallite.  It  is  evident  from  the  figure  that  the  linewidths  of  all  the  Raman 
modes  are  strongly  correlated  with  the  impurity  concentration.  The  values  of  the  Raman 
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linewidths  at  zero  slit  width  may  be  used  to  evaluate  the  phonons  lifetime,  t,  via  the  energy-time 
uncertainty  relation  [18], 

—  =  -  (2) 

h  t 

where  AE  is  the  Raman  linewidth  in  units  of  cm  *,  and  fi  =5.3*10  cm  -sec.  Table  I 
summarizes  the  results:  for  each  sample  the  zero  slit  linewidth  is  presented  along  with  the 
calculated  lifetime  for  each  of  the  Raman  modes.  Two  main  conclusions  can  be  drawn  from  the 

above  results:  the  lifetimes  are  correlated  with  impurity  concentration,  and  the  mode  of  both 

samples  has  a  significantly  longer  lifetime  than  those  of  the  E2,  El(TO),  Al(TO),  and  Al(LO) 

modes. 

The  first  finding  may  be  indicative  of  the  existence  of  a  phonon  lifetime  shortening 
mechanism  via  impurities  and  will  be  discussed  next.  One  fundamental  mechanism  which 
affects  phonon  lifetimes  in  semiconductors  is  the  anharmonic  decay  of  a  phonon  into  other 
Brillouin  zone  phonons  such  that  the  energy  and  momentum  is  conserved  in  the  process  [20]. 
However,  other  channels  of  phonon-decay  are  clearly  possible  especially  in  material  containing 
impurities  and  defects  [5,  21].  These  crystal  imperfections  can  destroy  the  translational 
symmetry  of  the  crystal,  and  thus  perturb  the  characteristic  lifetimes  of  the  propagating  phonons. 
Our  findings  indicate  that  both  channels  for  the  phonon  decay  exist  in  the  AIN  crystallites; 
therefore,  the  measured  lifetime,  x,  in  Equation  1  may  be  given  by  [18]. 
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(3) 


r  rA  r, 

where  ta  and  t|  are  the  characteristic  decay  times  due  to  the  crystal  anharmonicity  and  impurity 
scattering,  respectively.  It  is  not  a  straightforward  task  to  deconvolute  the  shortening 
mechanisms  of  the  phonon  lifetimes  but  the  following  approximate  treatment  gives  an  estimate 
for  ti.  Specifically,  due  to  its  high  crystal-quality  and  relatively  low  impurity  content,  the 
measured  phonon  lifetimes  of  the  transparent  crystallite  can  be  assumed  to  be  controlled  mainly 
by  the  anharmonic  interaction  and  thus  to  be  equal  to  the  value  of  tA  of  each  of  the  phonons. 
Using  the  values  of  tA  as  the  characteristic  anharmonic  decay  times  of  AIN,  the  decay  times  due 
to  impurities  in  the  blue  crystallite  can  be  calculated  via  Equation  3.  This  approximate  method 
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yields  values  of  xi  of  3.7  ps,  0.89  ps,  1.27  ps,  0.72  ps,  and  0.74  ps  for  the  E2  ,  E, ,  El(TO), 
Al(TO),  and  Al(LO)  modes,  respectively. 

The  anharmonic  decay  of  phonons  has  been  studied  extensively  for  zinc-blende 
semiconductors  [20,  29-30]  but  has  received  considerably  less  attention  for  wurtzite  structure 
materials  such  as  AlxGa,.xN  [23,  31].  Of  special  significance  for  the  nitride-based  wurtzite 
structures  is  the  decay  of  the  LO  phonon  modes  which  dominates  in  the  polar-optical-phonon 
carrier  interaction  in  many  high  speed  optoelectronic  devices  [32-33].  The  accepted  model, 
known  as  Kiemens’  model,  of  LO  phonon  decay  in  cubic  structure  materials  has  been  established 
to  occur  via  the  annihilation  of  the  LO  Raman  phonon  (of  frequency  coLO)  and  the  creation  of  two 
LA  phonons  (each  of  frequency  ©la)  at  the  zone  edge  [20].  In  this  model,  a  special  case  of  the 
three  phonon  process,  energy  conservation  requires  that  2©la=©lo-  However,  in  AIN,  due  to  the 
substantial  mass  difference  between  the  A1  and  the  N  elements  a  large  gap  in  the  AIN  phonon 
spectrum  exists:  the  Al(LO)  frequency  at  the  T  point  is  at  ~  890  cm'1  while  the  LA  at  zone  edge 
has  frequency  ~  400  cm'1  [34].  Accordingly  the  Klemens’  channel  is  not  applicable  here 
since  2©La  <  ©lo  •  A  four-phonon  interaction  would  permit  the  conservation  of  energy;  however, 
this  interaction  would  result  in  an  unreasonably  long  lifetime  on  the  order  of  two  or  three  orders 
of  magnitude  longer  than  that  predicted  for  the  three  phonon  process  [23].  To  account  for  the 
lifetime  of  the  LO  in  the  wurtzite  materials,  a  three-phonon  model  has  been  proposed  recently  in 
which  the  LO  mode  decays  into  a  TO  mode  and  an  LA  mode;  a  decay  route  which  results  in  a 
more  reasonable  lifetime  of  a  few  picoseconds  [23].  This  model  has  been  supported  by  time- 
resolved  Raman  measurements  of  the  Al(LO)  phonon  lifetime  in  wurtzite  GaN  thin  films  [31]. 
Moreover,  our  experimental  results  for  the  Al(LO)  lifetimes  of  the  AIN  crystallites  (~0.5  ps)  are 
in  accord  with  the  time  scale  of  the  three-phonon  anharmonic  process,  if  we  include  impurity 
scattering  as  an  additional  shortening  mechanism. 

In  order  to  further  investigate  the  relatively  long  lifetime  of  the  E\  mode,  the  Raman 
modes  of  GaN  as  well  as  ZnO  crystals  were  analyzed.  Figures  8  and  9  depict  the  Raman  spectra 
of  the  GaN  crystallite:  the  Al(TO),  El  (TO),  E2, ,  and  E2  modes  are  at  530  cm'1,  558  cm'1,  567 
cm'1,  and  142  cm'1,  respectively.  The  corresponding  lifetimes  computed  via  the  previously- 
described  method  are  summarized  in  Table  I.  Similar  to  the  AIN  lifetimes,  the  E2  mode  of  the 
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GaN  has  a  significantly  longer  lifetime  than  those  of  the  E*,  El(TO),  and  Al(TO)  modes. 
Figure  10  depicts  the  linewidth  behavior  of  the  e\  and  the  E;  modes  of  the  ZnO,  GaN,  and  the 
transparent  AIN  crystals;  the  zero  slit  lifetimes  of  the  ZnO  modes  are  presented  in  Table  I.  Our 
analyses  indicate  that  the  E,  mode  of  the  ZnO  crystallite  exhibits  a  relatively  long  lifetime  as 
well.  As  was  noted  previously,  the  LO  mode-lifetime  in  wurtzite  structure  has  been  considered 
[23].  However,  the  long  lifetime  of  the  El2  mode  in  that  structure  is  not  unexpected  in  view  of 
energy-conservation  constraints,  the  densities  of  the  final  states,  and  the  fact  that  the  anharmonic 
coefficient  has  not  been  determined  theoretically.  The  energy  of  the  E^  mode  lies  in  the  low 

energy  regime  of  the  wurtzite  dispersion  curves  [34-35]  and  only  the  acoustical  phonons  are 
available  as  the  final  states  of  the  decay  process.  Indeed  at  energies  where  the  densities-of-states 
are  appreciable  for  the  final  state  phonons,  the  energies  of  these  phonons  are  comparable  with 
those  of  the  Eo  mode;  thus,  it  may  be  reasonable  to  expect  that  the  matrix  elements  for  energy- 
conserving  decay  processes  multiplied  by  the  corresponding  densities-of-states  be  reduced, 
relative  to  the  other  modes.  Such  a  reduced  scattering  rate  would,  of  course,  result  in  a  longer 
lifetime  of  the  mode.  A  full  understanding  of  our  measured  value  for  the  lifetime  of  the  e\ 

mode  must  await  detailed  theoretical  calculations. 

Lastlv,  we  investigated  the  mode  of  GaN  film.  Fig.  1 1  depicts  the  Raman  spectra  of  a 

3  pm  undoped  GaN  film  grown  on  a  SiC  substrate  with  an  AIN  buffer  layer  and  with  a  free 
carrier  concentration  ~  10 16  cnr3  [7],  The  spectra  were  acquired  in  a  backscattering  geometry  at 
room  temperature  (RT)  and  at  10  K;  the  slit  width  was  set  at  200  pm.  The  lineshape  of  both  lines 
is  a  Lorentzian  with  a  linewidth  of  4  cm'1  and  3.2  cm  1  for  the  RT  and  the  10  K  spectrum, 
respectively.  Using  Equations  1  and  2,  the  zero  slit  values  and  the  lifetimes  may  be  calculated. 
The  RT  and  the  10  K  lifetimes  were  found  to  be  2  ps  and  1.5  ps,  respectively.  Our  results 
indicate  that  the  lifetime  of  the  E*  mode  of  the  GaN  film  is  similar  to  that  of  the  crystallite  (1.4 
ps).  Furthermore,  no  significant  lifetime  change  occurs  at  10  K.  A  weak  and  broad  Al(LO) 
Raman  line  from  the  film  was  observed  in  the  spectra:  however,  due  to  the  strong  SiC  Raman 
signal  which  appears  in  the  same  energy  range,  an  accurate  analysis  of  the  Al(LO)  lifetime  was 
not  possible.  The  lifetime  of  the  Al(LO)  mode  was  obtained  previously  [31]  via  time-resolved 
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Raman  measurements;  this  reported  lifetime  was  ~  3  ps.  This  lifetime  is  of  the  same  order-of- 
magnitude  as  the  lifetimes  reported  here  for  the  other  phonons  in  the  GaN  crystallite;  see  Table  I. 

IV.  CONCLUSIONS 

The  lifetimes  of  the  Raman  phonons  of  GaN  and  AIN  were  found  to  exhibit  two  time 
regimes:  the  long  lifetime  range  for  the  E-)  phonons  and  the  much  shorter  lifetime  range  for  the 
E' ,  El  (TO),  and  Al(TO)  phonons.  The  lifetime  of  the  Al(LO)  mode  of  AIN  crystallites  were 
found  to  be  the  shortest  among  the  modes  studied  here.  The  Ej>  lifetime  of  high-quality  and  of 
relatively  low-impurity-concentration  AIN  crystallite  is  4.4  ps,  which  is  about  five  times  longer 
than  that  of  the  El ,  El(TO),  and  Al(TO)  modes  and  about  an  order  of  magnitude  longer  then 
those  of  the  Al(LO)  mode.  The  lifetime  of  the  e\  mode  of  the  GaN  crystallite  is  10.1  ps,  which 

is  about  ten  times  longer  than  those  of  the  other  modes.  We  also  investigated  two  AIN 
crystallites  differing  in  their  impurity  concentration  and  found  that  the  phonon  lifetimes  were 
correlated  with  the  impurity  concentrations.  The  effect  of  the  two  order-of-magnitude  increase  in 
relative  impurity  concentration  shortened  the  lifetime  of  the  phonons  by  ~  50  /o.  Additionally, 
we  studied  the  E*  phonon  lifetime  of  a  high-quality  GaN  film  at  room  temperature  and  at  10  K. 
Our  results  indicate  that  the  phonon  lifetime  of  the  film  is  similar  to  that  of  the  crystallite, 
moreover,  no  significant  difference  between  the  lifetimes  at  the  two  temperature  ranges  was 

observed. 
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Table  I:  The  zero  slit-width  Raman  linewidths  and  the  corresponding  phonon  lifetimes  of  the 
AIN,  GaN,  and  ZnO  crystallites. 
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0.91 
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0.53 
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0.76 

14.4 
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0.46 
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11.9 

0.45 

18.9 

0.28 
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FIGURE  CAPTIONS 


Fig.  1.  A  micrograph  of  the  transparent  AIN  crystallite;  size  ~  500  x  50  pm. 

Fig.  2.  The  Raman  spectrum  (at  100  pm  slit  width)  of  the  Al(TO)  mode  of  the  AIN  crystallite. 
The  line  and  the  dashed  line  are  the  Lorentzian  and  the  Gaussian  fits  to  the  data,  respectively. 

Fig.  3.  The  Raman  spectra  and  the  curve  fit  to  the  spectra  of  the  transparent  AIN  crystallite 
(upper  figure).  The  deconvolved  spectra  is  presented  in  the  lower  figure. 

Fig.  4.  The  Raman  spectra  and  the  curve  fit  to  the  spectra  of  the  blue  AIN  crystallite  (upper 
figure).  The  deconvolved  spectra  is  presented  in  the  lower  figure. 

Fig.  5.  The  Raman  spectra  of  the  eJ>  mode  of  the  transparent  AIN  (lower  spectrum),  and  of  the 
blue  crystallite. 

Fig.  6.  The  Raman  spectra  of  the  Al(LO)  mode  of  the  transparent  AIN  (lower  spectrum)  and  of 
the  blue  crystallite. 

Fig.7.  The  linewidth  of  the  Raman  mode  as  a  function  of  the  slit  width.  The  squares  represent 
the  data  of  the  blue  AIN,  and  the  dots  of  the  transparent  AIN.  The  lines  are  the  fits  to  the  data 

using  Eq.  1 . 

Fig.  8.  The  Raman  spectra,  the  curve  fit,  and  the  deconvolution  spectra  of  the  GaN  crystallite. 

Fig.  9.  The  Raman  spectra  of  the  e\  mode  of  the  GaN  crystallite. 

Fig.10.  The  e!>  modes  (lower  curves)  and  the  E^  modes  (upper  curves)  of  the  blue  AIN 

(squares),  ZnO  (circles),  and  GaN  (diamonds)  as  a  function  of  the  slit  width. 

Fig.  1 1.  The  spectra  (200  pm  slit  width)  of  E‘  modes  of  GaN  film  at  room  temperature  and  at 
10  K.  The  dashed  line  represents  the  Lorentzian  fits  to  the  data. 
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ABSTRACT 

Raman  analysis  of  the  E2  mode  of  AlxGai.xN  in  the  composition  range  0  <  x  <  1  is 
presented.  The  lineshape  was  observed  to  exhibit  a  significant  asymmetry  and  broadening  toward 
the  high  energy  range.  The  spatial  correlation  model  is  discussed,  and  is  shown  to  account  for  the 
lineshape.  The  model  calculations  also  indicate  the  lack  of  a  long-range  order  in  the  CVD 
(chemical  vapor  deposition)  alloys.  These  results  were  confirmed  by  X-ray  scattering:  the  relative 
intensity  of  the  superlattice  line  was  found  to  be  negligible.  The  line  broadening  of  the  E2  mode 
was  found  to  exhibit  a  maximum  at  a  composition  x~0.5  indicative  of  a  random  disordered  alloy 
system.  The  stress  state  of  the  alloys  was  found  to  be  tensile  and  was  attributed  to  the  difference  in 
the  thermal  expansion  coefficients  of  the  SiC  substrate  and  the  film. 

INTRODUCTION 

Raman  spectroscopy  has  been  proven  in  the  past  to  be  a  powerful  method  in  studying  alloy 
systems  among  which  are  AlxGai.xAs,  InxGai.xAs,  and  many  more  [1-2].  Raman  scattering  from 
a  perfect  crystal,  which  consists  of  a  lattice  with  translational  symmetry,  involves  phonons  with 
wavevector  q=0  near  the  Brillouin  zone  (BZ)  center  [3].  The  alloying  may  destroy  the 
configurational  translational  symmetry,  and  the  phonons  of  q>0  may  be  allowed  to  participate  in 
the  Raman  scattering.  As  a  result  the  Raman  lineshape  changes,  and  an  investigation  of  the  line 
characteristics  may  yield  valuable  information  on  the  microstructure  state  of  the  alloy. 

In  this  paper  we  present  observations  and  analytical  study  of  an  asymmetric  E2  Raman 
lineshape  of  AlxGai.xN  films.  The  asymmetry  of  the  line  is  shown  to  agree  with  the  spatial 
correlation  model  which  is  based  on  the  relaxation  of  the  q=0  selection  rules  due  to  the  alloying. 
The  calculations  indicated  that  the  AlxGai_xN  lattice  lacks  a  long  range  order,  a  result  which  was 
confirmed  by  X-ray  diffraction:  the  superlattice  diffraction  line  was  found  to  exhibit  a  very  weak 
relative  intensity.  Moreover,  the  functional  behavior  of  the  E2  linewidth  with  respect  to  the 
composition  x  in  the  range  0-1  indicated  a  random  distribution  of  the  alloy  constituents  with 
maximum  disorder  at  composition  x~0.5. 

In  general,  stress  may  be  analyzed  via  the  Raman  peak  position.  However,  the  peak 
position  of  the  alloy  depends  on  the  convoluted  effects  of  the  stress,  the  relaxation  of  the  q=0 
selection  rules,  and  the  mode-type.  These  effects  are  functions  of  the  composition;  the  stress  is 
also  a  function  of  the  lattice  constants  and  the  thermal  expansion  coefficients  of  the  film-substrate 
system.  Thus,  in  a  sample  of  known  spatially  uniform  composition,  and  with  structural  defects 
due  to  stress,  the  stress-type  may  be  analyzed  by  acquiring  spectra  at  a  defect  and  at  a  defect-free 
region.  Our  analysis  using  this  approach  indicates  that  the  AlGaN  alloys  are  under  tensile  stress. 
The  stress  is  suggested  to  be  due  to  the  difference  in  thermal  expansion  coefficients  of  the  film  and 
the  SiC  substrate. 

EXPERIMENT 

Raman  spectra  were  acquired  in  a  backscattering  geometry  in  which  the  E2  is  an  allowed 
Raman  mode  from  the  c-face  of  the  Wurtzite  structure  [4],  In  the  experiments  the  488  nm  line  of 
an  Ar+  ion  laser  was  utilized  and  the  data  were  acquired  at  sample  temperatures  of  300K  and  10K. 
The  data  for  the  lineshape  investigation  were  acquired  in  a  macro  mode  of  laser  spot  size  ~  2  mm; 


for  the  stress  measurements  the  micro  mode  of  spot  size  ~  3  (im  was  utilized.  The  experimental 
error  in  the  linewidth  and  peak  position  measurements  is  ±  1.5  cm*1  and  0.5  cm*1  respectively. 
The  AlxGai_xN  films  were  grown  via  the  organometallic  chemical  vapor  deposition  (OMCVD) 
method  at  ~  1100  °C  on  6H-SiC(0001)  substrates  with  a  100  nm  AIN  buffer  layer.  The  thickness 
of  the  films  is  ~  2  |im  and  the  composition,  x,  which  was  determined  via  Rutherford  backscattering 
(RBS),  energy  dispersive  X-ray  (EDX),  and  Auger  spectroscopy,  is:  0.06,  0.12,  0.22,  0.32,  and 
0.70.  The  Raman  data  points  for  x=0  and  x=l  were  obtained  from  GaN  film  and  AIN  crystallite 
respectively. 

RESULTS  AND  DISCUSSION 

The  Configurational  Disorder  in  AlGaN  Films 

Figure  1  shows  the  room-temperature  Raman  spectra  of  the  E2  line  from  AlxGai.xN  films 
of  compositions  0.06,  0.12,  0.22,  and  0.70,  which  exhibit  linewidths  T  of  8,  13,  16,  and  19  cm*1 
respectively.  As  depicted  in  the  figure,  the  spectral  lineshape  for  films  of  x>0.06  exhibit 
asymmetric  broadening  and  a  peak  shift  toward  higher  frequency.  Possible  line  broadening 
mechanisms  applicable  to  alloys  include  thermal  broadening,  activation  of  a  symmetry  forbidden 
zone-center  (q=0)  mode  which  lies  in  the  same  frequency  range  as  the  investigated  line,  and 
broadening  due  to  activation  of  a  collection  of  modes  of  wavevectors  q>0.  The  last  two 
broadening  mechanisms  result  from  the  elimination  of  the  translational  symmetry  of  the  lattice  due 
to  alloying. 

To  investigate  the  thermal  contribution  to  the  line  broadening,  Raman  spectra  were  acquired 
at  T=10K,  and  no  significant  change  was  observed  between  the  shapes  and  linewidths  of  the  room 
and  cold  temperature  spectra.  Thus  the  linewidth  in  our  sample  is  not  strongly  affected  by 
temperature.  The  only  effect  of  the  low  temperature  is  the  shifting  of  the  peak  position  by  ~  2  cm*1 
toward  higher  frequency,  which  was  also  previously  observed  in  GaN  films  [5]  and  crystals  [6] 
and  was  attributed  there  to  the  thermal  contraction  of  the  bonding. 

The  second  possible  mechanism  is  the  activation  of  a  q=0  symmetry  forbidden  mode  which 
might  be  convoluted  with  the  E2  line  and  cause  the  asymmetric  broadening.  However,  the  only 
mode  in  the  frequency  range  of  the  E2  line  that  is  forbidden  in  the  back-scattering  geometry  is  the 
Al(TO)  mode  of  GaN  at  -560  cm*1  which  lies  at  a  lower  frequency. 

A  more  plausible  mechanism  to  account  for  the  high  frequency  asymmetric  linewidth  in  our 
Raman  spectra  is  the  spatial  correlation  model,  also  known  as  the  confinement  model.  The  model 
was  developed  to  explain  the  asymmetry  line  in  BN  [7]  and  Si  [8]  and  has  been  successfully 
applied  to  quantify  the  lineshape  behavior  of  Gai_xAlxAs  and  Gaj.xInxAs  alloy  systems  [1].  The 
foundation  of  the  model  lies  in  attributing  the  relaxation  of  the  q=0  Raman  selection  rules  to  the 
phonon  confinement  in  a  finite  domain  of  size  L.  The  size  of  L  in  an  alloy  system  may  be  viewed 
as  the  average  size  of  the  ordered  domains  which  are  embedded  in  the  configurational-disordered 
matrix.  According  to  the  model,  as  L  gets  smaller,  the  range  of  wavevectors  Aq  becomes  larger:  a 
wider  range  of  frequencies  are  allowed  to  participate  in  the  Raman  scattering,  and  as  a  result  the 
Raman  lineshape  exhibits  a  change.  The  lineshape-change  depends  on  the  phonon  dispersion 
curve  co(q);  for  a  curve  of  negative  slope  the  line  asymmetry  is  toward  lower  frequencies  and  for  a 
positive  slope  the  asymmetry  is  toward  higher  frequencies.  The  peak  position  also  shifts 
accordingly.  However,  the  position  is  a  function  of  the  stress  and  the  alloy-mode  as  well:  an  alloy 
can  be  categorized  as  being  a  one-  or  two-  mode  type  [9]  and  the  behavior  of  the  peak  position 
depends  on  the  type.  The  deconvolution  of  the  peak  position  is  left  to  a  future  investigation;  at 
present  we  address  merely  the  origin  of  the  asymmetric  lineshape  and  line  broadening. 

The  Raman  intensity  I(co)  for  the  spatial  correlation  model  may  be  expressed  as  [1, 7-8] 


where  q  is  in  units  of  2jt/a,  a  is  the  lattice  constant,  and  To  is  the  linewidth  of  the  material  of 
composition  x=0  which  is  GaN  in  this  study.  Our  GaN  films  exhibit  linewidths  in  the  range  of 


3-5  cm'1  depending  on  the  film  quality.  For  the  calculations  we  used  the  upper  range  value  to 
account  for  the  stress  in  the  alloy  which  was  inferred  from  the  appearance  of  cracks  in  the  films. 
In  the  calculation  we  assumed  a  spherical  correlation  domain  and  BZ. 

To  evaluate  Equation  1  we  fitted  the  phonon  dispersion  curves,  which  have  been  calculated 
by  Azuhata  et  al  [10],  in  the  vicinity  of  BZ  center.  In  their  calculations  the  E2  mode  along  the  A 
direction  exhibits  a  doubly  degenerate  behavior  and  is  a  weak  function  of  q,  thus  having  no  effect 
on  the  line  broadening  mechanism.  In  contrast,  the  E2  curves  along  the  I  and  T  directions  both 
split  into  a  lower  and  an  upper  branch,  where  the  former  is  a  decreasing  function  of  q  and  the  later 
is  an  increasing  function;  these  curves  are  expected  to  influence  the  Raman  lineshape.  Our 
calculations  of  Eq.  1  best  concur  with  the  phonon  dispersion  curve  of  the  form:  co(q)  =  A  +  Bq2 
(in  units  of  cm"1  and  with  A=568  and  B=100)  which  is  the  approximated  behavior  of  the  upper 
branch  of  the  I  and  T  curves  in  the  vicinity  of  the  zone  center.  However,  when  we  included  the 
contribution  to  the  lineshape  arising  from  the  lower  branches  of  the  dispersion  curves  the  fit  was 
not  ideal.  Similar  behavior  has  been  previously  observed  in  the  calculation  of  BN  lineshape  [7]. 
The  reason  for  the  dominance  affect  of  the  upper  branches  in  our  calculations  is  under 
investigation. 

The  calculated  lines  and  the  Raman  data  for  samples  of  composition  0.12  and  0.22  are 
presented  in  Figure  1 .  As  depicted  in  the  figure  the  asymmetric  Raman  data  concur  with  the  spatial 

correlation  model,  where  the  small  deviation  ~  560  cm-1  may  be  attributed  to  the  activated  Al(LO) 
mode.  The  ordered  domain  size,  L,  was  determined  from  our  calculations  to  be  ~  2.5  -  3.5  nm;  a 
relatively  small  value  suggesting  the  lack  of  a  long  range  order  in  the  alloy  system.  In  order  to 
further  investigate  this  hypothesis  a  film  of  composition  x=0.22  was  characterized  by  XRD:  the 
scan  is  presented  in  Figure  2.  The  figure  shows  the  diffraction  peak  of  the  (0002)  planes  of  the 
AlGaN  (superimposed  is  also  the  SiC  substrate  peak)  and  the  inset  presents  the  weak  (0001) 
diffraction  peak.  The  (0001)  is  a  forbidden  reflection  in  the  WZ  structure  and  appears  as  a 
superlattice  line  only  when  lattice  ordering  exists  [11];  the  intensity  ratio  (0001)/(0002)  is  a 
measure  of  the  extent  of  the  order  [12].  AlxGai-xN  films  which  were  grown  by  molecular  beam 
epitaxy  have  been  reported  by  Korakakis  et  al  [1 1]  to  exhibit  a  long  range  order;  the  intensity  ratio 
of  their  XRD  lines  is  -  0.05  in  contrast  to  our  much  smaller  intensity  ratio  of  -  0.002.  Thus  the 
XRD  confirms  the  Raman  analysis,  indicating  that  the  alloys  in  this  study  do  not  exhibit  significant 
long  range  ordering.  More  insight  into  the  alloy  micro-structure  may  be  obtained  from  the  Raman 
linewidth  behavior.  Figure  3  presents  the  Raman  linewidth  of  the  E2  mode  as  a  function  of  the 
composition.  The  linewidth  behavior  follows  a  pattern  which  exhibits  a  maximum  at  composition 
x~0.5,  a  value  at  which  a  maximum  disorder  should  be  expected  in  a  random  system.  Thus  our 
data  suggest  that  the  distribution  of  the  alloy  constituents  for  the  most  part  is  random. 

The  Stress  State  of  the  AlGaN  Films 

The  morphology  of  our  AlGaN  films  consists  of  uniformly  distributed  crack  density  which 
is  indicative  of  the  high  stress  state  present  in  the  films.  In  general,  stress  may  be  analyzed  via  the 
shifting  of  the  Raman  peak  position.  However,  as  mentioned  the  Raman  peak  position  of  the  alloy 
depends  on  the  combined  effects  of  the  stress,  the  relaxation  of  the  qsO  selection  rules,  and  the 
mode-type.  All  these  effects  are  functions  of  the  composition;  the  stress,  however,  is  also  a 
function  of  the  lattice  constants  and  the  thermal  expansion  coefficients  of  the  film-substrate  system. 

In  order  to  determine  the  stress-type  in  the  alloy  we  took  advantage  of  the  film  morphology: 
Micro-Raman  spectra  were  acquired  at  a  crack  region,  where  the  film  is  relatively  relaxed,  and  at  a 
region  remote  from  the  crack,  where  the  film  is  strained.  Figure  4  shows  characteristic  spectra  of 
the  relaxed  and  strained  regions.  The  linewidth  of  both  spectra  is  —16  cm-1  implying  that  no 
significant  spatial  variation  of  composition,  x,  occurs  across  the  film.  Since  the  composition  is 
similar  in  the  relaxed  and  the  strained  regions,  the  Raman  peak  position  is  influenced  mainly  by  the 
stress  effects  due  to  the  film-substrate  thermal  and  lattice  mismatch.  As  can  be  observed  in  the 
figure,  the  peak  position  of  the  strained  region  exhibits  a  shift  of  ~3  cm-1  towards  lower 
frequencies  relative  to  the  peak  of  the  relaxed  region  suggesting  that  the  original  stress  state  of  the 
film  (before  relaxation  via  cracking  occurs)  is  tensile. 

A  source  of  the  tensile  stress  in  the  alloys  may  be  the  difference  in  the  thermal  expansion 
coefficients,  a,  of  the  film,  the  buffer-layer,  and  the  substrate  at  the  elevated  temperature  (1400  K) 


at  which  our  films  were  grown.  The  a’s  of  GaN  and  SiC  at  this  temperature  range  have  not  yet 
been  measured.  However,  at  -750  K  the  thermal  expansion  coefficients  along  the  a  axis  (aa)  and 
c  axis  (etc)  have  been  measured  [13];  for  GaN,  aa=6.2  and  ac=6.1,  and  for  SiC  aa=4.2  and 
ac=4.0,  where  a  is  in  units  of  10*6  K'1.  The  average  thermal  expansion  coefficients  of  AIN  at 
750  K  and  1200  K  have  been  reported  to  be  5.8  and  6.2  10‘6  K'1  respectively  [14].  The  smaller 
value  of  the  SiC  thermal  coefficient  relative  to  that  of  the  alloy  constituents  thus  may  be  the  cause 
of  the  tensile  stress  in  our  films.  To  acquire  further  insight  into  this  issue,  the  Raman  spectra  of 
AIN  film,  grown  under  similar  conditions  as  the  alloys,  was  compared  to  that  of  AIN  crystallite. 
The  spectra  are  presented  in  Figure  5  (the  spectrum  of  the  crystallite  was  acquired  with  a 
backscattering  geometry  from  the  a-face  for  which  the  El(TO)  is  an  allowed  mode  as  well).  As 
can  be  seen  in  the  figure,  the  AIN  Raman  peak  of  the  film  exhibits  a  tensile  shift  of  -4  cm-1  relative 
to  the  peak  position  656  cm-1  of  the  crystal.  A  result  suggestive  of  the  presence  of  stress  in  the 
film  due  to  the  different  coefficients  of  thermal  expansion.. 

CONCLUSIONS 

The  asymmetric  behavior  of  the  E2  Raman  lineshape  of  AlxGai-xN  alloys  was  attributed  to 
the  activation  of  phonons  of  q>0  arising  from  the  disordered  state  of  the  alloys.  The  line 
broadening  of  the  E2  Raman  mode  exhibits  a  maximum  at  composition  x~0.5  which  is  consistent 
with  a  random  distributed  system.  X-ray  diffraction  supports  these  findings:  the  superlattice 
relative  line-intensity  is  weak.  The  alloy  exhibits  tensile  stress  which  is  suggested  to  be  due  to 
interfacial  thermal  mismatch. 
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FIG.  1 .  The  Raman  spectra  of  the  E2  mode  of  films  of  various  composition  x. 
The  solid  lines  were  calculated  from  the  spatial  correlation  model  (Eq.  1). 


FIG.  2.  The  (0002)  x-ray  diffraction  of  film  of  x=0.22.  FIG.  3.  The  E2  Raman  linewidth  as  a  function 
The  inset  shows  the  very  weak  (0001)  superlattice  line,  of  composition  x. 


Raman  Intensity  (A.U.) 


540  550  560  570  580  590  600  610 

Raman  Shift  (cm' *) 


FIG.  4.  The  Raman  spectra  of  the  E2  mode  acquired  at  a  relaxed  and  at  a  strained  region 
of  AlGaN  film  of  composition  x=0.22. 
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FIG.  5.  Raman  spectra  of  (a)  AIN  film  grown  on  SiC  substrate  and  (b)  AIN  crystal. 
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ABSTRACT 

Raman  analysis  of  the  E2  mode  of  AlxGai.xN  in  the  composition  range  0  <  x  <  1  is  presented. 
The  lineshape  was  observed  to  exhibit  a  significant  asymmetry  and  broadening  toward  the  high 
energy  range.  The  spatial  correlation  model  is  discussed,  and  is  shown  to  account  for  the 
lineshape.  The  model  calculations  also  indicate  the  lack  of  a  long-range  order  in  the  chemical  vapor 
deposition  (CVD)  alloys.  These  results  were  confirmed  by  X-ray  scattering:  the  relative  intensity  of 
the  superlattice  line  was  found  to  be  negligible.  The  line  broadening  of  the  E2  mode  was  found  to 
exhibit  a  maximum  at  a  composition  x  =  0.5  indicative  of  a  random  disordered  alloy  system. 
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The  realization  of  achieving  high  quality  AlxGai_xN  films  has  prompted  intensive 
investigation  concerning  device  fabrication  and  material  and  device  properties.  Raman 
spectroscopy  has  been  proven  in  the  past  to  be  a  powerful  method  in  studying  alloy  systems 
among  which  are  AlxGa].xAs.  InxGai_xAs,  CaxSr].xF2,  ZnSi-x  Sex,and  many  more  [1-5], 
Raman  scattering  from  a  perfect  crystal,  which  consists  of  a  lattice  with  translational  symmetry, 
involves  phonons  with  wavevector  q^O  near  the  Brillouin  zone  (BZ)  center  [6].  The  alloying  may 
destroy  the  configurational  translational  symmetry,  and  the  phonons  of  q>0  may  be  allowed  to 
participate  in  the  Raman  scattering.  As  a  result  the  Raman  lineshape  changes,  and  an  investigation 
of  the  line  characteristics  may  yield  valuable  information  on  the  microstructure  state  of  the  alloy. 

In  this  Letter  we  present  for  the  first  time  observations  and  analytical  study  of  an 
asymmetric  E2  Raman  lineshape  of  AlxGai-xN  films.  The  asymmetry  of  the  line  is  shown  to  agree 
with  the  spatial  correlation  model  which  is  based  on  the  relaxation  of  the  qsO  selection  rules  due  to 
the  alloying.  The  calculations  indicated  that  the  AlxGai.xN  lattice  lacks  a  long  range  order,  a  result 
which  was  confirmed  by  X-ray  diffraction  (XRD):  the  superlattice  diffraction  line  was  found  to 
exhibit  a  very  weak  relative  intensity.  Moreover,  the  functional  behavior  of  the  E2  linewidth  with 
respect  to  the  composition  x  in  the  range  0  -  1  indicated  a  random  distribution  of  the  alloy 
constituents  with  maximum  disorder  at  composition  x=0.5. 

Raman  spectra  were  acquired  in  a  backscattering  geometry  in  which  the  E2  and  the  Al(LO) 
are  the  allowed  Raman  modes  from  the  c-face  of  the  Wurtzite  (WZ)  structure  [7].  Unlike  the  non¬ 
polar  E2  mode,  the  Al(LO)  mode  possesses  an  electric  field  which  interacts  with  that  produced  by 
the  free  carrier  waves  known  as  plasmons.  For  samples  with  a  sufficiently  high  free  carrier 
concentration,  this  phonon-plasmon  interaction  affects  the  Al(LO)  Raman  lineshape  significantly 
and  results  in  a  Raman  signal  with  very  weak  intensity  [8].  The  free  carrier  concentration  in  the 
samples  studied  here,  as  measured  via  the  Hall  effect,  is  relatively  high  (on  the  order  of 
5xl017/cm3)  and  as  a  result  the  Al(LO)  mode  is  not  detectable  in  the  Raman  spectra.  The  lack  of 
the  electric  field  interaction  makes  the  E2  mode  ideal  for  the  present  study.  Future  investigations 
may  address  the  lineshape  of  the  Al(LO)  mode. 
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In  the  experiments  the  488  nm  line  of  an  Ar+  ion  laser  was  utilized  and  the  data  were 

acquired  at  sample  temperatures  of  300K  and  at  10K.  The  experimental  error  in  the  linewidth  and 

peak  position  measurements  is  ±  1.5  cm'1  and  0.5  cm'1  respectively.  The  AlxGai_xN  films  were 

grown  via  the  organometallic  chemical  vapor  deposition  (OMCVD)  method  at  ~  1 100  °C  on  6H- 

o 

SiC(OOOl)  substrates  with  a  1000  A  AIN  buffer  layer  [9].  The  thickness  of  the  films  is  ~  2  [im 
and  the  composition,  x,  which  was  determined  via  Rutherford  backscattering  (RBS),  energy 
dispersive  X-ray  (EDX),  and  Auger  spectroscopy,  is:  0.06,  0.12,  0.22,  0.32,  and  0.70.  The  data 
points  for  x=0  and  x=l  were  obtained  from  GaN  film  and  AIN  crystallite  [10]  respectively. 

Figure  1  shows  the  room-temperature  Raman  spectra  of  the  E2  line  from  AlxGai_xN  films 
of  compositions  0.06,  0.12,  0.22,  and  0.70,  which  exhibit  linewidths  of  8,  13,  16,  and  19  cm'1 
respectively.  As  depicted  in  the  figure,  the  spectral  lineshape  for  films  of  x>0.06  exhibit 
asymmetric  broadening  and  a  peak  shift  toward  higher  frequency.  Possible  line  broadening 
mechanisms  applicable  to  alloys  include  thermal  broadening,  activation  of  a  symmetry  forbidden 
zone-center  (q=0)  mode  which  lies  in  the  same  frequency  range  as  the  investigated  line,  and 
broadening  due  to  activation  of  a  collective  of  modes  of  wavevectors  q>0.  The  last  two 
broadening  mechanisms  result  from  the  elimination  of  the  translational  symmetry  of  the  lattice  due 
to  alloying. 

To  investigate  the  thermal  contribution  to  the  line  broadening,  Raman  spectra  were  acquired 
at  T=10K,  and  no  significant  change  was  observed  between  the  shapes  and  linewidths  of  the  room 
and  cold  temperature  spectra.  Thus  the  linewidth  in  our  sample  is  not  strongly  affected  by 
temperature.  The  only  effect  of  the  low  temperature  is  the  shifting  of  the  peak  position  by  ~  2  cnr1 
toward  the  high  frequency,  which  was  also  previously  observed  in  GaN  films  [8]  and  crystals  [11] 
and  was  attributed  there  to  the  thermal  contraction  of  the  bonding. 

The  second  possible  mechanism  is  the  activation  of  a  q=0  symmetry  forbidden  mode  which 
might  be  convoluted  with  the  E2  line  and  cause  the  asymmetric  broadening.  However,  the  only 
mode  in  the  frequency  range  of  the  E2  line  that  is  forbidden  in  the  back-scattering  geometry  is  the 
Al(TO)  mode  of  GaN  at  -560  cm'1  [7]  which  lies  at  a  lower  frequency. 
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A  more  plausible  mechanism  to  account  for  the  high  frequency  asymmetric  linewidth  in  our 
Raman  spectra  is  the  spatial  correlation  model,  also  known  as  the  confinement  model.  The  model 
was  developed  to  explain  the  asymmetry  line  in  BN  [12]  and  Si  [13]  and  has  been  successfully 
applied  to  quantify  the  lineshape  behavior  of  Gai-xAlxAs  and  Gai_xInxAs  alloy  systems  [1].  The 
foundation  of  the  model  lies  in  the  wavevector  uncertainty:  Aq=2rt/L  where  Aq  is  the  phonon 
wavevector  range,  L  is  a  spatial  dimension  into  which  the  phonon  is  confined,  and  q  is  related  to 
the  Raman  frequency  co  via  the  phonon  dispersion  relation  to(q).  The  size  of  L  in  an  alloy  system 
may  be  viewed  as  the  average  size  of  the  ordered  domains  which  are  embedded  in  the 
configurational-disordered  matrix.  According  to  the  model,  as  L  gets  smaller  Aq  becomes  larger:  a 
wider  range  of  frequencies  are  allowed  to  participate  in  the  Raman  scattering,  and  as  a  result  the 
Raman  lineshape  exhibits  a  change.  The  lineshape-change  depends  on  the  phonon  dispersion 
curve  0)(q);  for  a  curve  of  negative  slope  the  line  asymmetry  is  toward  lower  frequencies  and  for  a 
positive  slope  the  asymmetry  is  toward  higher  frequencies.  The  peak  position  also  shifts 
accordingly.  However,  the  position  is  a  function  of  the  alloy-mode  as  well:  an  alloy  can  be 
categorized  as  being  a  one-  or  two-  mode  type  [14-15]  and  the  behavior  of  the  peak  position 
depends  on  the  type.  The  deconvolution  of  the  peak  position  is  left  to  a  future  investigation;  at 
present  we  address  merely  the  origin  of  the  asymmetric  lineshape  and  line  broadening. 

The  Raman  intensity  I(co)  for  the  spatial  correlation  model  may  be  expressed  as  [1, 12-13] 


where  q  is  in  units  of  2rt/a,  a  is  the  lattice  constant,  and  To  is  the  linewidth  of  the  material  of 
composition  x=0  which  is  GaN  in  this  study.  Our  GaN  films  exhibit  linewidth  in  the  range  3  -  5 
cm'1  depending  on  the  film  quality  [8].  For  the  calculations  we  used  the  upper  range  value  to 
account  for  the  stress  in  the  alloy  which  was  inferred  from  the  appearance  of  cracks  in  the  films. 
In  the  calculation  we  assumed  a  spherical  correlation  domain  and  BZ. 
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To  evaluate  Equation  1  we  fitted  the  phonon  dispersion  curves,  which  have  been  calculated 
by  Azuhata  et  al  [16],  in  the  vicinity  of  BZ  center.  In  their  calculations  the  E2  mode  along  the  A 
direction  exhibits  a  doubly  degenerate  behavior  and  is  a  weak  function  of  q,  thus  having  no  effect 
on  the  line  broadening  mechanism.  In  contrast,  the  E2  curves  along  the  X  and  T  directions  both 
split  into  a  lower  and  an  upper  branch,  where  the  former  is  a  decreasing  function  of  q  and  the  later 
is  an  increasing  function;  these  curves  are  expected  to  influence  the  Raman  lineshape.  Our 
calculations  of  Eq.  1  best  concur  with  the  phonon  dispersion  curve  of  the  form:  co(q)  =  A  +  Bq2 
(in  units  of  cm'1  and  with  A=568  and  B=100)  which  is  the  approximated  behavior  of  the  upper 
branch  of  the  X  and  T  curves  in  the  vicinity  of  the  zone  center.  However,  when  we  included  the 
contribution  to  the  lineshape  arising  from  the  lower  branches  of  the  dispersion  curves  the  fit  was 
not  ideal.  Similar  behavior  has  been  previously  observed  in  the  calculation  of  BN  lineshape  [12]. 
The  reason  for  the  dominance  affect  of  the  upper  branches  in  our  calculations  is  under 
investigation. 

The  calculated  lines  and  the  Raman  data  for  samples  of  composition  0.12  and  0.22  are 

presented  in  Figure  1 .  As  depicted  in  the  figure  the  asymmetric  Raman  data  concur  with  the  spatial 

correlation  model,  where  the  small  deviation  ~  560  cm'1  may  be  attributed  to  the  activated  Al(LO) 

o 

mode.  The  ordered  domain  size,  L,  was  determined  from  our  calculations  to  be  ~  25  -  35  A ;  a 
relatively  small  value  suggesting  the  lack  of  a  long  range  order  in  the  alloy  system.  In  order  to 
further  investigate  this  hypothesis  a  film  of  composition  x=0.22  was  characterized  by  XRD:  the 
scan  is  presented  in  Figure  2.  The  figure  shows  the  diffraction  peak  of  the  (0002)  planes  of  the 
AlGaN  (superimposed  is  also  the  SiC  substrate  peak)  and  the  inset  presents  the  weak  (0001) 
diffraction  peak.  The  (0001)  is  a  forbidden  reflection  in  the  WZ  structure  and  appears  as  a 
superlattice  line  only  when  lattice  ordering  exists  [17];  the  intensity  ratio  (0001)/(0002)  is  a 
measure  of  the  extent  of  the  order  [18].  AlxGai.xN  films  which  were  grown  by  molecular  beam 
epitaxy  have  been  reported  by  Korakakis  et  al  [17]  to  exhibit  a  long  range  order;  the  intensity  ratio 
of  their  XRD  lines  is  ~  0.05  in  contrast  to  our  much  smaller  intensity  ratio  of  ~  0.002.  Thus  the 
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XRD  confirms  the  Raman  analysis,  indicating  that  the  alloys  in  this  study  do  not  exhibit  significant 
long  range  ordering. 

More  insight  into  the  alloy  micro-structure  may  be  obtained  from  the  Raman  linewidth 
behavior.  Figure  3  presents  the  Raman  linewidth  of  the  E2  mode  as  a  function  of  the  composition. 
The  linewidth  behavior  follows  a  pattern  which  exhibits  a  maximum  at  composition  x=0.5,  a  value 
at  which  a  maximum  disorder  should  be  expected  in  a  random  system.  The  solid  line  in  the  figure 
was  calculated  (in  arbitrary  units)  using  the  relation  for  the  entropy  of  mixing  which  is  a  measure 
of  the  cofigurational  disorder  in  a  random  alloy  [19].  Thus  our  data  suggest  that  the  distribution  of 
the  alloy  constituents  for  the  most  part  is  random.  The  issue  of  randomness  in  AxBi-xC  alloys 
has  been  previously  addressed  and  it  was  suggested  that  the  anions,  C,  are  periodically  distributed 
in  their  sub-lattice  sites  while  the  cations,  A  and  B,  are  randomly  distributed  in  the  other  sub-lattice 
[2,  20].  However,  achieving  an  ordered  alloy  in  some  families  of  ternary  tetrahedral 
semiconductors  has  been  demonstrated  to  be  a  function  of  the  growth  temperature  [20];  the 
conditions  for  the  ordered  state  in  the  AlxGai-xN  system  have  yet  to  be  explored. 

To  conclude,  the  asymmetric  behavior  of  the  E2  Raman  linewidth  of  AlxGai.xN  alloys  was 
attributed  to  the  activation  of  phonons  of  q>0  arising  from  the  disordered  state  of  the  alloys.  X-ray 
diffraction  supports  the  finding;  the  supperlattice  relative  line-intensity  is  weak.  The  line 
broadening  of  the  E2  Raman  mode  exhibits  a  maximum  at  composition  x=0.5  which  is  consistent 

with  a  random  distributed  system. 
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FIGURE  CAPTIONS 

Figure  1.  The  RT  Raman  spectra  of  the  E2  mode  of  AlxGai-xN  films  of  composition  x=0.06, 
0.12, 0.22  and  0.70.  The  respective  linewidths  are:  8, 13,  16,  and  19  cnr1,  and  peak  positions  at: 
567,  569,  572,  579  cm'1.  The  solid  lines  superimposed  on  the  spectra  of  films:  x=0.12  and  0.22 
were  calculated  from  the  spatial  correlation  model  presented  in  Eq.  1. 

Figure  2.  The  (0002)  X-ray  diffraction  of  Alo.22Gao.7sN;  the  scan  displays  the  SiC  substrate 
peak.  The  inset  to  the  figure  shows  the  very  weak  (0001)  superlattice  line. 

Figure  3.  The  E2  Raman  linewidth  as  a  function  of  composition  (dots).  The  solid  line  is  obtained 
from  calculations  of  the  entropy  of  mixing  of  an  alloy  system  (in  arbitrary  units). 
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